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Abstract. We review recent experimental and theoretical work on superconductivity in 
ultrasmall metallic grains, i.e. grains sufficiently small that the conduction electron energy 
spectrum becomes discrete. The discrete excitation spectrum of an individual grain can be 
measured by the technique of single-electron tunneling spectroscopy, and reveals parity effects 
indicative of pairing correlations in the grain. After introducing the discrete BCS model 
that has been used to model such grains, we review a phenomenological, grand-canonical, 
variational BCS theory describing the paramagnetic breakdown of these pairing correlations 
with increasing magnetic field. We also review recent canonical theories that have been 
developed to describe how pairing correlations change during the crossover, with decreasing 
grain size, from the bulk limit to the limit of few electrons, and compare their results to those 
obtained using Richardson's exact solution of the discrete BCS model. 
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1 Introduction 

Since its discovery by Kammcrlingh Chines in 1911, superconductivity has become 
one of the most-studied phenomena in condensed matter physics, and its microscopic 
explanation via the highly successful pairing theory proposed in 1957 by Bardeen, 
Cooper and Schrieffer (BCS) [@ is one of the landmark achievements of 20th century 
physics. Yet, despite the long history of the subject, to this day experimental advances 
in sample fabrication and measurement techniques continue to reveal novel aspects 
of superconductivity, which often require extensions or modifications of the existing 
theoretical framework. 

The subject of this review, superconductivity in ultrasmall metallic grains, is a 
case in point: In the mid 1990's, Ralph, Black and Tinkham (RBT) succeeded for the 
first time to directly measure the discrete excitation spectrum of individual ultrasmall 
metallic grains (of radii r < 5 nm and mean level spacings d > 0.1 meV) using 
a technique called single-electron-tunneling spectroscopy: by attaching such a grain 
via oxide tunnel barriers to two leads they constructed a single-electron transistor 
having the grain as central island, and showed that a well-resolved, discrete excitation 
spectrum could indeed be extracted from the conductance 0-0. 

This opened up a new frontier in the study of electron correlations in metals, since 
the ability to resolve discrete energy levels allows the nature of electron correlations 
to be studied in unprecedented detail. Since 1995, RBT's technique has been used 
to probe superconducting pairing correlations in Al grains fi [?|, noncquilibrium ex- 
citations 0-[ffo| and spin-orbit interactions p, [|, [ll], |l2|, [l3| in normal grains, and 
ferromagnetic correlations in Co grains llS]. A comprehensive survey of all ex- 
perimental and theoretical developments (up to March 2000) relating to spectroscopic 
studies of ultrasmall metallic grains may be found in the review of von Delft and 
Ralph |]l6) . The present review is an excerpt of Ref. Jll| , and is devoted exclusively to 
superconductivity (more precisely, to pairing correlations) in ultrasmall metallic grains. 

For several reasons, RBT's experiments on superconducting pairing correlations in 
ultrasmall Al grains attracted quite some attention (T^-ETJ: 

First, for largish > 5 nm) Al grains, RBT's measurements revealed a rather 
striking parity effect |j 0, [l?]] : a grain with an even number of electrons had a distinct 
spectroscopic gap d) but an odd grain did not. This is clear evidence for the 
presence of superconducting pairing correlations in these grains, and indicates that a 
BCS-like theory would be appropriate for their description. 

Second, the spectroscopic gap for even grains was driven to zero by an applied 
magnetic field, hence the paramagnetic breakdown of pairing correlations could be 
studied in detail ||, |§. 

Third, even the "largish" grains were so small that standard grand-canonical BCS 
mean field theory is no longer applicable: (a) the single-particle mean level spacing 
d = 1/A/"(ef) ~ 1/Vol [where Af(e) is the density of states per spin species] was 
comparable to the bulk superconducting gap, which we shall denote by A, so that a 
mean field approach is no longer reliable (it requires d -C A); and (b), the number of 
electrons on such a grain is well-defined, hence a canonical theory is required. RBT's 
experiments stimulated the development of corresponding extensions of BCS theory. 

Fourth, in RBT's smallest grains (r < 3 nm), the distinct spectroscopic gap ob- 
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served for largish even grains could no longer be unambiguously discerned. This ob- 
servation revived an old but fundamental question: What is the lower size limit for the 
existence of superconductivity in small grains?: Anderson had addressed this question 
already in 1959 |4l|, arguing that if the sample is so small that its electronic eigenspec- 
trum becomes discrete, "superconductivity would no longer be possible" when its mean 
level spacing d becomes larger than the bulk gap A. Heuristically, this is obvious (see 
Fig. [7] below) : d is the number of free-electron states that pair-correlate (those with 
energies within A of ep), i.e. the "number of Cooper pairs" in the system; when this 
becomes < 1, it clearly no longer makes sense to call the system "superconducting". 

Although Anderson's answer is correct in general, it generates further questions: 
What, precisely, does "superconductivity" mean in ultrasmall grains, for which many 
of the standard criteria such as zero resistivity, Meissner effect and Josephson effect, 
are not relevant^? What happens in the regime d > A in which superconductivity has 
broken down? Is the breakdown parity dependent? How is it influenced by a magnetic 
field? This review attempts to provide detailed answers to these and related questions. 

The review is divided into two distinct parts. Part I (Sees. || to ^|) is devoted 
to RBT's experiments and their detailed theoretical interpretation. After briefly dis- 
cussing the experimental setup and summarizing the main experimental results, we 
analyze and qualitatively explain the latter in the framework of a phenomenologi- 
cal theory by Braun et al. [[l8| [p| ^2|. This theory offers a simple intuitive picture 
for visualizing the pairing correlations and how these change when the grain size is 
decreased. 

Part II (Sees. [l(] to |TJ) is devoted to further theoretical developments, inspired 
by RBT's experiments but not directly concerned with their interpretation Jl7|-p0[. 
In particular, RBT's new experiments stimulated a number of theoretical attempts to 
quantitatively describe the crossover from the bulk limit d« A, where superconduc- 
tivity is well-developed, to the fluctuation-dominated regime of d S> A, where pairing- 
correlations survive only in the form of weak fluctuations. Describing this crossover 
constituted a conceptual challenge, since the standard grand-canonical mean-field BCS 
treatment of pairing correlations Jl7]]-[^4| breaks down for d > A. This challenge 
elicited a series of increasingly sophisticated canonical treatments of pairing correla- 
tions J29J -J37j, based on a simple reduced BCS Hamiltonian for discrete energy levels, 
which showed that the crossover is completely smooth, but, interestingly, depends on 
the parity of the number of electrons on the grain |]l7| . Very recently, the main con- 
clusions of these works were confirmed |35) using an exact solution of the discrete-level 
BCS model, discovered by Richardson in the context of nuclear physics in the 1960s 
[^3)-(5^]. (The existence of this solution came as a surprise - in the form of a polite 
letter from its inventor - to those involved with ultrasmall grains, since hitherto it had 
apparently completely escaped the attention of the condensed- matter community.) 

A detailed outline of the contents of the two parts may be found in the opening 
paragraph of each, or in the table of contents. 



a For an isolated nm-scale grain, (i) its resistivity is not defined, since electron motion is ballistic 
and the mean free path is boundary-limited; (ii) the grain radius is smaller than the penetration 
depth, so that no Meissner effect occurs; and (iii) the electron number is fixed, so that the order 
parameter cannot have a well-defined phase. 
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Part I: Experiment and phenomenological theory 

Part I of this review is devoted to RBT's experiments and their detailed theoretical 
interpretation. It is organized as follows: section 

(2) is devoted to experimental details; 

(3) summarizes RBT's main experimental results; 

(4) proposes a phenomenological model for an isolated ultrasmall grain; 

(5) discusses how pairing correlations can be visualized in a fixed- N system and 
explains when and in what sense it can be called "superconducting" ; 

(^) presents a generalized variational BCS approach for calculating the eigenenergies 

of various variational eigenstates of general spin s; 
(0) discusses how an increasing magnetic field induces a transition from a pair- 
correlated state to a normal paramagnetic state; 
(||) presents theoretical tunneling spectra of the RBT type, which are in qualitative 
agreement with RBT's measurements; 

discusses various parity effects that are expected to occur in ultrasmall grains. 



2 Experimental details 

In the first generation of RBT's experiments of 1995 j^, ^, Q], a grain made from Al 
(a superconducting material) was connected to two metal leads via high-resistance 
tunnel junctions, with capacitances C L and C R , say. In the next generation of 1997 
[0, the grain was also coupled capacatively to a gate, with capacitance C g . The 
resulting device, schematically depicted in Fig. |l|(a), has the structure of a SET, with 
the grain as central island. The circuit diagram for an SET is shown in Fig. |l|(b). 
Applying a bias voltage V between the two leads causes a tunnel current / to flow 
between the leads through the grain, via incoherent sequential tunneling through the 
tunnel junctions. The current can be influenced by changing the gate voltage V g 
(hence the name "transistor"), which tunes the electrostatic potential on the grain 
and thereby also its average number of electrons N. (For devices without a gate these 
two quantities cannot be tuned and instead have some sample-dependent, fixed value. 
For such devices, set C g = in all formulas below.) 




G 



Ann. Phys. (Leipzig) 10 (2001) 3 



The physics of SETs had been clarified in the early 1990s |52| through extensive 
studies of lithographically defined SETs of mesoscopic size, i.e. with micron-scale 
central islands. The fundamentally new aspect of RBT's work was that their SETs, 
made by a novel fabrication technique (described in Ref. B), were nanoscopic in size: 
they had ultrasmall grains with radii between 15nm and 2nm as central islands, which 
were thus several orders of magnitude smaller in volume than in previous experiments. 
This had two important consequences: 

1. The grain's charging energy Eq = e 2 /2C was much larger than for mesoscopic 
SETs, ranging roughly between 5 and 50 meV (where C = C L + C R + C g ). Eq is 
the scale that determines the energy cost for changing the electron number N by one. 
Since for ultrasmall grains it far exceeds all other typical energy scales of the SET, 
such as those set by the bias voltage (V < 1 mV), the temperature (T < 4.2 K) 
and the bulk superconducting gap for Al (Abuik = 0.18 meV), fluctuations in electron 
number are strongly suppressed. 

2. Discrete eigenstates of the conduction electron energy spectrum became resolvable 
- their mean level spacing d ranged from 0.02 to 0.3 meV. This agrees in order of 
magnitude with the estimate d = 1/A/"(ef) obtained using the free-electron expression 
for the density of states A/"(£f) at the Fermi surface of a 3D grain, 

2n 2 h 2 _ 1.50eV.nm 2 
~ m/cpVol ~~ k F Vol ' ^ ' 

where, for example, k-p = 17.5 nm" 1 for Al. The measured d-values are much larger 
than k^T for the lowest temperatures attained (around T ~ 30mK), but on the order 
of Abuik- However, the number of conduction electrons for grains of this size is still 
rather large (between 10 4 and 10 5 ). 

Since the two scales Eq and d differ by at least an order of magnitude, they 
manifest themselves in two distinct and easily separable ways in the low-temperature 
I-V curves of RBT's devices: 

1. Coulomb-blockade phenomena: When V is varied on a large scale of tens of 
mV for fixed V g [Fig. §, the I-V curves have a typical "Coulomb-staircase form" 
characteristic of SETs: zero current at low \V\ (the "Coulomb blockade" regime), 
sloping steps equally spaced in V, and step thresholds sensitive to V g . This proves 
that the tunnel current flows only through one grain. The maximal width of the flat 
step of zero current around \eV\ = is governed, in order of magnitude, by Ec and 
typically varies between 5 and 50 mV. 

As V^, is varied, the I-V curves periodically repeat. The "orthodox theory" for 
Coulomb blockade phenomena ]54, 55, 56|, [l(| explains this as follows: the electrostatic 



work required to add iV ox excess electrons with a total charge of eA^x (with e < 0) 
to a grain with initial random off-set charge Q , while the time- independent voltages 
Vl, and V g of the left and right leads and the gate electrode, respectively, are held 
fixed, has the form 

£ P ot(iVex) = eV N ex + E c Nl . (2) 
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Fig. 2 Current-voltage curves for an ultrasmall SET at 50 mK, artificially offset on the 
vertical axis for a set of equally-spaced values of V g between —1.2 and 1.8 V. The I-V curves 
display Coulomb-staircase structure on a bias voltage scale of tens of mV. By fitting these 
to standard SET theory p^|, the SET's basic parameters can be determined: = 3.5 Mf2, 
R R = 0.2 MQ,, C L = 3.5 aF, C R = 9.4 aF, C g = 0.09 aF, E c = 46 meV. Assuming the 
grain shape to be hemispherical and using the measured capacitances to estimate its surface 
areas (see J?], |l6| for details), the grain's radius and mean level spacing were estimated as 
r ~ 4.5 nm and d ~ 0.45 meV. 



Here Vd = (Qo + H r =L,R,g Cr^J /C represents the electrostatic potential on the 

grain, and EqN^ x represents the Coulomb interaction energy of the N ex excess elec- 
trons due to their mutual repulsion. Since, for given V g , the system adjusts N ex such 
as to minimize E pot (N ex ), the I-V characteristics are V^-periodic, with period e/Cg. 

2. Fine- structure due to level discreteness: When V is varied on the much smaller 
scale of a few mV around the threshold of the Coulomb blockade regime and the 
temperature is sufficiently low (T <C d), the I-V curves have a step-like substructure, 
and correspondingly, the differential conductance (dl/dV) curves contain a series of 
fine peaks, see Fig. |[ As first pointed out by Averin and Korotkov |57| , such small 
steps in the I-V curve are expected to occur whenever the voltage drop across one 
of the tunnel junctions (say r = L, R) equals the threshold energy at which the rate 
for tunneling across that junction into or out of one of the grain's discrete energy 
eigenstates becomes nonzero, since this opens up another channel for carrying current 
across that junction. Such tunneling processes are ilustrated schematically in Fig. 0. 



More formally, let S^* Q , denote the rate for the tunneling transition \a'] 



N' 



\o)n between two grain eigenstates (with electron number N' and N = N' ± 1, and 
eigenenergies , £^ , respectively,), induced by transferring an electron across barrier 
r onto (upper sign) or off (lower sign) the grain. A golden rule calculation (see Sec. 2.3.3 
of Rcf. |l6| for details) then shows that these rates have the form 

= fK - C T e(V r - V D ))r r ± a , , (3) 

where f(E) = l/(e E ^ kBT + 1) is the Fermi function, V r — V^> is the voltage drop 
(electrostatic potential difference) between lead r and the grain, and are transition 
probabilities. With increasing transport voltage V = V^ — V&, a current step thus 
occurs each time a "bottleneck rate" (associated with a rate-limiting tunneling 
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Fig. 3 Excitation spectra for the same sample J?j] as in Fig. |2| measured at T = 50 mK 
and H — 0.05 T (to drive the Al leads normal), for four different V g -values, corresponding to 
different values for the grain's average electron number (from top to bottom: N + l, N, N, N — 
1). The curves are artificially offset on the vertical axis and each is labeled by the associated 
bottleneck tunneling rate (from initial state |G)jv' to final state |a)jv, two of which 

are illustrated schematically in Fig. ^), the bottleneck barrier being r = L in this case. 
Plotted is dl/dV vs. energy, where the latter is given by the voltage-to-energy conversion 
factor |eV|(CR + C g /2)/C = 0.73|eV|, which reflects the voltage drop across barrier L (for 
a derivation of this factor, see Sec. 2.3 of Ref. |^|). The sizeable spectroscopic gap between 
the first two peaks in the middle two curves, and its absence in the top and bottom curves, 
reflects the pairbreaking energy cost in the excitation spectrum of a superconducting grain 
with an even number of electrons, and implies that N is odd. 



process) is switched from "off" (exponentially small) to "on" (of order T^,), i.e. each 
time e(V r — Vrj) passes through a threshold at which one of the inequalities 

±e(V r -V )> (4) 

becomes true. The tunneling spectrum yields particularly useful information if the 
initial state for all bottleneck tunneling processes is always the ground state, \a') n' = 
\G)n' (whose electron number N' is determined by V g ). This will be the case if the 
following conditions are met: (i) for the given (fixed) value of gate voltage V g , all 
bottleneck processes involve the same^] barrier; (ii) the temperature is sufficiently low 
(T <C d); and (iii) nonequilibrium effects are negligible (requiring relaxation rates on 
the grain to be much greater than tunneling rates). Under these conditions, which are 
satisfied by the data of Fig. || (and Fig. || below), the distances between the current 
steps or conductance peaks directly reflect the so-called fixed-N excitation spectrum 
of the grain, i.e. the set of energy differences 

between the eigenenergies £^ of all those iV-electron eigenstates \o)n that are acces- 
sible final states from the initial ground state |G)jv<, via a tunneling processes onto 
the grain if JV = N' + 1, or off the grain if N = N' - 1. 

For devices having a gate, two very interesting options exist: Firstly, by tuning V g 

b If bottleneck rates for tunneling across both barriers are comparable, the probability of finding N' 
and JV electrons on the grain will be comparable, and the analysis is considerably more complicated; 
see Sec. 2.3.4 of Jlil for details. 
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N is odd 




Fig. 4 Schematic depiction of bottleneck tunneling processes governing the excitation spec- 
tra of Fig. for TV being odd, with rates: (a) t G N ana - 00 ^a^G N+ i' corresponding 
to two different choices of V s just below or above the degeneracy point at which the TV- and 
(TV + l)-electron ground states are degenerate. The long-dashed line indicates the equilib- 
rium, V = chemical potential of the L and R leads. Solid (dashed) arrows depict bottleneck 
tunneling transitions into the lowest- (highest) energy final states accessible for the chosen 
value of V , and filled circles represent the electron configuration of the lowest-energy final 
state. (For a more detailed discussion of such diagrams, in particular how they change with 
V g , see Sec. 2 of Ref. @.) 



such that the Coulomb blockade regime is large or small, so that the ^-threshold at 
which current begins to flow is large or small, nonequilibrium effects can be maximized 
or minimized, respectively, depending on whether one chooses to study them or not. 
In this review, we consider only the near-equilibrium situation. Strong nonequilibrium 
was studied in Refs. 0-Jl(| and is reviewed in Sec. 6 of Ref. Secondly, by tuning 
V g by an amount large enough (~ E^/e) to change TV' by one unit, the influence on 
the spectrum of the parity of the number of electrons on the grain can be studied. 

Parity effects occur, for instance, in the magnetic-field dependence of the fixed- 
TV excitation spectrum, which can be obtained by simply tracing the motion of the 
conductance peak positions as a magnetic field is turned on (at fixed V^). This is 
shown in Fig. below, which nicely illustrates the breaking of Kramers degeneracy 
by the applied field: For H = 0, the grain will have time-reversal symmetry. For an 
even-TV grain, the many-electron wave function for the ground state will be a spin 
singlet, in order that the orbital energy be minimized. In contrast, the ground state 
of an odd- TV grain for H = necessarily is two-fold degenerate, by Kramers' theorem, 
forming a Kramers doublet. When H is turned on, this doublet is Zeeman-split by 
±^/iB<7-ff. Therefore, for an even-N grain at small H, the lowest-lying tunneling 
excitations correspond to transitions from the even-grain ground state singlet to the 
odd-grain ground state doublet, i.e. to two states split by H, so that the lowest-^ 
conductance peak will exhibit Zeeman splitting in an applied field [Fig. [1(b)] . On the 
other hand, for an odd-N grain with T <C /iB 9 H/kB, the odd-grain ground state will 
be the lower-energy state of the Kramers doublet; the lowest-lying tunneling excitation 
will thus consist only of a single transition from this odd-grain ground state to the 
even-grain ground state singlet, so that the lowest-]/ conductance peak will not split 
into two as a function of H [Fig. 0(a)]. It follows that in Fig. ||, TV is odd. 
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Fig. 5 Magnetic field dependence Q of excitation spectra such as those of Fig. ^| and taken 
for the same grain, at (a) 14 ~ HO mV and (b) V g ~ 180 mV [corresponding to the second 
and topmost curves of Fig. B, respectively] . Each line represents a distinct conductance peak 
in the dl/dV curve and traces how its energy changes with H. Upward- moving peaks are 
broader and less distinct than downward-moving ones (for reasons poorly understood) and 
can be followed only for a limited range of H before they are lost in the background. The 
distances between lines directly give the grain's (a) fixed-(7V + 1) and (b) fixed-iV excita- 
tion spectrum, where N is odd; the corresponding bottleneck processes are schematically 
illustrated in Figs. ^](a) and (b), respectively. The vertical dashed lines indicate the first 
four level-crossing fields H s<s , (defined in Eq. @ and assigned by comparison with Fig. [l^, 
see Sec. |), namely H ,i = 4T, #1/2,3/2 = 4.25T, Hi, 2 = 5.25T and #3/2,5/2 = 6.5T with 
uncertainty ±0.13T (half the ^-resolution of 0.25T). 



3 A gap in the excitation spectrum 

The spectra shown in Figs. || and [| which are typical for RBT's results for largish 
Al grains (r > 5 nm), reveal a very striking feature: if the number of electrons on 
the grain in the final state after the bottleneck tunneling process is even (middle two 
curves of Fig. |^), the excitation spectra display a spectroscopic gap between the first 
two conductance peaks that is significantly larger than the mean spacing between 
subsequent peaks, whereas no such gap occurs if the final-state electron number is 
odd (top and bottom curves of Fig. ||). In other words, even- TV excitation spectra 
(number parity p = 0) are gapped, whereas odd-TV excitation spectra (p = 1) are 
not. This is even more clearly apparent when comparing Figs. ^|(a) and ^(b), which 
give the magnetic-field (H) dependence of an even-A~ and odd-A~ excitation spectrum, 
respectively. However, in their smallest grains (r < 3 nm) no such clear parity- 
dependent spectroscopic gap could be discerned. 

BRT interpreted these observations as evidence for superconducting pairing corre- 
lations in their larger grains, using notions from the BCS theory of superconductivity: 
in an even grain, all excited states involve at least one broken Cooper pair, i.e. two 
BCS quasiparticles, and hence lie at least 2A above the fully-paired BCS ground 
state; in contrast, in an odd grain all states have at least one unpaired electron, i.e. 
at least one quasiparticle, and hence no significant gap exists between ground- and 
excited states. Fig. [|is a cartoon illustration of this interpretation, if one imagines two 
electrons drawn on the same energy level to represent a "Cooper pair" (making this 
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cartoon precise will be one of the main goals of review): in Fig. [|(a) the final electron 
number is even, and all final excited states (reached via dashed arrows) have one less 
"Cooper pair" than the final ground state (reached via the solid arrow); in contrast, 
in Figs. 0(b) the final electron number is odd, and the final ground and excited states 
have the same number of "Cooper pairs" . 

The approximately linear ff-dependence of the excitation spectra in Fig. |5| was 
attributed by RBT to the Zeeman energy shifts of discrete levels in a magnetic field (its 
effect on orbital motion is neglible, as will be discussed in Sec. fh^ ). The fact that the 
lowest state in Figs. ||(a) or (b) does not or does display Zeeman splitting, respectively, 
implies that N is odd, as mentioned above. The reduction of the spectroscopic gap 
in Fig. ||(a) therefore is purely due to Zeeman energy shifts and has nothing to do 
with the reduction of the BCS gap parameter due to pair-breaking that occurs in bulk 
samples in a magnetic field [psj . A detailed discussion of the spectra's magnetic field 
dependence will be given in Sees. ^ and |^. 

For completeness, it should be remarked that a spectral gap in ultrasmall super- 
conducting grains was observed as long ago as 1968 by Giaever and Zeller J59J, po[ , 
who studied tunneling through granular thin films containing electrically insulated Sn 
grains. They found gaps for grain sizes right down to the critical size estimated by 
Anderson (radii of 2.5 nm in this case), but were unable to prove that smaller parti- 
cles are always "normal" . RBT's experiments are similar in spirit to this pioneering 
work, but their ability to focus on individual grains makes a much more detailed study 
possible. 



4 A discrete BCS model for ultrasmall grains with pairing correlations 

In this section we construct a model for an isolated ultrasmall grain with pairing 
correlations, using phenomcnological arguments valid for the regime d < A. The 
model, which we shall call "discrete BCS model" , allows a rather detailed qualitative 



understanding of the measurements of RBT 18, 19| and hence is "phenomcnologically 



successful". For d ^> A it is unrealistically simple, however, and should rather be 
viewed as a toy model for studying how pairing correlations change as the grain is 
made smaller and smaller. 



4.1 A simple reduced BCS interaction plus a Zeeman term 

Following the philosophy of the "orthodox theory" for Coulomb blockade phenomena 
[ p4| . |l6| , we assume that the only effect of the Coulomb interaction is to contribute an 
amount EqN^ [cf. Eq. (||)] to the eigenenergy of each eigenstate of the grain. Since 
the charging energy is huge (5 to 50 meV) in ultrasmall grains, this term strongly sup- 
presses number fluctuations around the optimal value of N ex , so that, to an excellent 
approximation, all energy eigenstates will also be number eigenstates.| Since EqN^ x 
is constant within every fixed- N Hilbert space, we shall henceforth ignore it, with the 
understanding that the model we are about to construct should always be solved in 

c An exception occurs at a so-called degeneracy point, where E po t(N ex ) = E po t(N ex + 1); the 
resulting complications will not be considered here. 
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a fixed- N Hilbert space (and that the use of grand-canonical approaches below, after 
having dropped EqN^ x , simply serves as a first approximation to the desired canonical 
solution) . 

The only symmetry expected to hold in realistic, irregularly-shaped ultrasmall 
grains at zero magnetic held is time-reversal symmetry. We therefore adopt a single- 
particle basis of pairs of time-reversed states |j±), enumerated by a discrete index 
j. Their discrete energies £j are assumed to already incorporate the effects of im- 
purity scattering and the average of electron-electron interactions, etc. As simplest 
conceivable model that incorporates pairing interactions and a Zeeman coupling to a 
magnetic held, we adopt a Hamiltonian H = Hq + H rc d of the following reduced BCS 
form: 

^0= ( £ J ah ) C )a C ]a . #rcd = ~Ad^ 4+4- C j- C j+ ■ ( 6 ) 

j,<r=± ij 

Here — ah = a^iiggH is the Zeeman energy of a spin a electron in a magnetic held 
H, and we shall take h > below. Models of this kind had previously been studied 
by Strongin et al fH, Miihlschlegel et al. ||, || and Kawataba g, ||. The hrst 
application to RBT's grains for h — was by von Delft et al. |T^] and for h ^ by 
Braun et al. @,|lJJ[l|. 

Due to level repulsion the Sj's will, to hrst approximation, be uniformly spaced. 
Unless otherwise specified, we shall for simplicity always (except in Sec. |l3]) take a 
completely uniform spectrum with level spacing d. Fluctuations in the level spacings 
have been studied with methods of random matrix theory p0|, with qualitatively 



similar results (see Sec. 13). For a system with a total of N electrons, where the 
electron number parity p = A r mod2 is equal to for even N and 1 for odd N, we 
use the label j = for the lowest-lying non-doubly-occupied level (with occupation 
number p) in the T = Fermi sea, which we shall denote by \Fn). We choose the 
Fermi energy at £f = writo^] 

s j =jd+(l-p)d/2, (7) 

thereby taking the doubly-occupied and empty levels of |Fjv) to lie symmetrically above 
and below £f (see e.g. Fig. || below)— The parameter [i in Eq. (||) is, in g.c. theories, 
the chemical potential, whose valueH determines the average particle number. For 
canonical theories, which make no reference to a chemical potential, fi is not needed 
and can be dropped (i.e. set equal to 0). 

The pairing interaction is of the reduced BCS form, in that it scatters a pair of 
electrons from one pair of time-reversed states into another .[] It is taken to include 
only states whose energy separation from the Fermi energy lies within the cutoff given 
by the Debye frequency: \sj\ < cjd- The pair-coupling constant in Eq. (^|) is written as 

d This convention differs slightly from that used in Jr7[ |f8j [hJ, namely Ej = jd + eq. The latter 
is a little less convenient, resulting, e.g., in a p-dependent chemical potential for the variational BCS 
ground states discussed below, /ip CS = so + (p — l)d/2, whereas (^) results simply in /ip CS = 0. 

°Note that the use of a reduced BCS interaction means that couplings between non-time-reversed 
pairs of states are neglected. A theoretical motivation, based oil random matrix theory, for this 
reduced form may be found in Ref. |10[ , or in Sec. 6.1.3 of Ref. Jig . Experimental evidence for the 
sufficiency of the reduced form is discussed in Refs. jl8| |l9| and, in most detail, in Sec. 4.7 of Ref. Eq|. 
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Xd, where A is a dimensionless parameter independent of the grain's volume, to make 
it explicit that both Hq and i? re d make extensive (cx Vol) contributions to the ground 
state energy (since the number of terms in each sum ^ . in Eq. (^|) scales with N, and 
d cx £f/N). The "bulk gap" of the model, obtained by solving the standard BCS gap 
equation [Eq. J27|)] at T = in the bulk limit, thus is 

A = w D /sinh(l/A) . (8) 

To be precise, by "bulk limit" we shall always mean d/A — ► and N — > oo while the 
product Nd is kept fixed, and use dJ2j ~ * f ^ £ j- 

An applied magnetic field will completely penetrate an ultrasmall grain, since its 
radius (typically r < 5nm) is much smaller than the penetration length of 50 nm for 
bulk Al. The Zeeman term in Eq. (^|) models the fact that RBT's measured tunnel 
spectra of Fig. H evolve approximately linearly as a function of magnetic field, with g 
factors between]] 1.95 and 2 (determined from the differences between measured slopes 
of up- and downward-moving lines). Deviations from g = 2 probably result from spin- 
orbit scattering, known to be small but nonzero in thin Al films |36|, [37| , but neglected 
below (where g = 2 is used). 

Intuitively speaking, it is clear that the discrete BCS model introduced above 
contains all ingredients necessary to make contact with the spectra of Fig. ||: it is 
formulated in terms of discrete levels, it contains a pairing interaction which is known, 
from bulk BCS theory 13, Bs| , to cause a gap in the excitation spectrum, and it contains 
a Zeeman term that will cause eigenenergies to linearly depend on an applied magnetic 
field. Indeed, we shall see in Sec. || that it can be used to obtain a rather detailed 
qualitative understanding of the spectra of Fig. |[ 



4.2 Why orbital diamagnetism is negligible in ultrasmall grains 

Of course, a magnetic field in principle also couples to the orbital motion of the 
conduction electrons - in bulk samples, this is the origin of the Meissner effect. Orbital 
effects in spherical and cylindrical superconductors whose dimensions are smaller than 
the penetration depth were first considered by Larkin f38| . However, in grains as small 
as those of RBT, orbital diagmagnetic effects are negligible plf , just as for thin films 
in a parallel magnetic field |6^, |6^|. The reason is as follows: 

Let -fforb denote the field scale above which orbital diamagnetism becomes impor- 
tant. If, in a random-matrix description of the grain's spectrum, iJ or b is formally 
associated with the field at which the crossover between the symplectic and unitary 
ensembles, driven by the orbital effects of the magnetic field, is complete, it is found 
@ that 

#orb » - = , (9) 

H N /.E Thouloss /d 



f Claims of smaller g factors made in Ref. [y| are wrong, the resu t of confusing different orbital 
states as Zeeman-split spin states. This was made clear in Ref. V\, where it was observed that 
upward-trending Zeeman states can have significantly smaller amplitude than downward-trending 
states, making them difficult to observe. 
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where Erhouiess is the Thouless energy and <E>o (= hc/2e — 2067T.nm 2 ) is the flux 
quantum. An intuitive understanding for the origin of this result can be obtained 
by the following argument (70): Associate H or b with the field at which the orbital 
splitting of the eigenenergies of two time-reversed states |j±) becomes comparable 
to the mean level spacing. Denoting the angular momentum of these states by 
±(l z )jK, the orbital diamagnetic contribution to their eigenenergies is ±(l z }j/jiBH, 
hence if rb w d/ {2{l z ) j^l-q). Now, the angular momentum of an electron travers- 
ing the "closed trajectory" corresponding to a discrete quantum level can be esti- 
mated as (l z )h ~ rn(Atypd/h), where the bracketed factor is the typical (directed) 
area At yp covered by its trajectory divided by the period h/d of its motion. The 
number of bounces off the grain's boundaries during this time is roughly i/rhouiess/rf, 
since the ft/i/rhouicss is the time to cross the grain once. Hence the directed area is 
A typ « r 2 \J E'YYiouiess/d, where the square root accounts for the fact that the direction 
of motion after each bounce is random [7l| ] . Collecting the various estimates results 
in Herb ~ $o / i^r 2 -&rhouicss/cO j which, up to a factor of 7r, agrees with Eq. (^). 

Now, assuming ballistic electron motion in the grain, the Thouless energy has the 
form i^Thouicss ~ hvF j (air) , where a is a geometrical constant of order unity. Using d 
from Eq. ([!]), we see that -fforb grows like r~ 3 with decreasing grain size. Taking a = 3 
(as in Ref. iQ), we find from Eq. (||) that hemispherical Al grains with radii of (say) 
r w 3 or 5 nm have iJ or b ~ 19 or 7 T, respectively. If larger values are used for a, as 
would be appropriate for more pancake-shaped grains ||, i? rb would be even larger. 

We may thus conclude that orbital diagmagnetic effects only begin to play a role 
for largish grains (> 5 nm), and then only for the highest fields (of 7 T) studied by 
RBT. Indeed, some larger grains do show slight deviations from H-linearity || for 
large fields, which probably reflect the onset of such orbital effects; however, these 
are much smaller than Zccman effects in the grains of present interest, and will be 
neglected here. Thus, the discrete BCS model assumes that Pauli paramagnetism due 
to the Zeeman energy completely dominates orbital diamagnetism, similarly to the 
case of thin films in parallel magnetic fields |66|, |67j. 

4.3 Choice of numerical values for model parameters 

When doing numerical calculations for the discrete BCS model, some choices must 
be made for the numerical values of its parameters (though slight changes in their 
values will not change the results qualitatively). We shall follow the choices made by 
Braun et al. | fl9| , since these led to reasonable agreement between experimental and 
theoretical excitation spectra. For the Debye frequency they used the textbook value 
[fr2[ for Al of lot) — 34meV. Making an appropriate choice for the "bulk gap" A is less 
straightforward, since its experimental value for systems of reduced dimensionality 
often differs from that of a truly bulk system, presumably due to (poorly-understood) 
changes in the phonon spectrum and the effective electron-phonon coupling. For ex- 
ample, for thin Al films f3l| , f73[ it is known that A t hinfiim — 0.38 meV, which is about 
twice as large as the gap of a truly bulk system, Abuik = 0.18 meV. (This increase in 
A is not universal, though; e.g., for Nb A is smaller in thin films than in the bulk.) 
Since ultrasmall grains are in many ways analogous to thin films in a parallel mag- 
netic field [see Sec. Q], Braun et al. adopted the thin-film value for grains too, i.e. 
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used A ~ 0.38meV. These choices imply that the dimensionless pair-coupling con- 
stant A = [smlT^o/A)]- 1 [cf.Eq. (§)] has the value A = 0.194. (In Sec. | we shall 
see, a posteriori, that the choices A = 0.34 and A = 0.189 would have been slightly 
more appropriate.) Finally, for those numerical calculations that are explicitly cut-off 
dependent, Braun et al. smeared the cutoff of the BCS interaction over two single- 
electron levels; this smooths out small discontinuities that would otherwise occur in 
d-dependent quantities each time the energy \ej\ of some large-|j| level moves beyond 
the cutoff lu-£) when d is increased. 

Note that the above way of choosing A lumps into a single phenomenological con- 
stant all the poorly-understood effects of reduced dimensionality |nj on the phonons 
that mediate the attractive electron-electron interaction. Studying these effects in de- 
tail would be interesting in its own right, but would require systematic investigations 
with grains of well-controlled shapes and sizes. For the case of RBT's irregularly- 
shaped grains, using a phenomenological coupling constant seems the best one can do. 
Note, though, that the precise value of A is not very important as long as all ener- 
gies are measured in units of A (as we shall do for all numerical calculations), since 
most of the A-dependence is thereby normalized away. Therefore, the slight difference 
between the A-values proposed above and those used in ]3l], [l2| ^5) (namely 0.224) 
hardly matters. 



4.4 Some general properties of the eigenstates - the blocking effect 

The eigenstates of the discrete BCS model of Eq. have some simple but general 
properties that are worth stating at the outset. 

Firstly, every eigenstate of H will also be an eigenstate of the number operator 

Secondly, since the interaction only involves levels within the cutoff energy cod of 
£f, the dynamics of those lying outside this range is trivial. We shall thus ignore them 
henceforth and focus only on the remaining set of interacting levels, denoting this set 
by 7. 

Thirdly, singly- occupied levels do not participate in the pair scattering described 
by H: "unpaired" electrons in such levels are not scattered to other levels, hence the 
labels of singly-occupied levels are good quantum numbers. Moreover, every unpaired 
electron Pauli-blocks the scattering of other pairs into its own singly-occupied level, 
i.e. it restricts the phase space available to pair scattering and thereby weakens the 
amount of pairing correlations, as we shall see in detail later. This was called the 
"blocking effect" by Soloviev j74|, who discussed it extensively in the early 1960's in 
the context of nuclear physics. The eigenstates |a) and corresponding eigenenergies 
£ a of H thus have the following general forms: 

\a) = |*„,S) = n c tj^>, (10) 

U n 
l*n) = £ V(jl,---,Jn)n & ]j VaC > - 
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Fig. 6 A cartoon depiction of the exact ground state for a reduced BCS Hamiltonian, for 
N being even (a) or odd (b): they are coherent superpositions of eigenstates of Ho (whose 
respective amplitudes are not depicted) that all have the same N; the leftmost of these is 
(a) the even or (b) the odd Fermi ground state |Fjv), whose Fermi energy is indicated by the 
wavy line. 



£ a = £ n +£B(h), £ B (h) =^2(e t - n- o-th) . (12) 

This describes N = 2n + b electrons, b of which are unpaired and sit in a set B of 
singly-occupied, blocked levels, making a contribution £s(h) to the total eigenenergy. 
The remaining n pairs of electrons, created by the pair operators b^ — ct_|_ej_, are 
distributed among the remaining set U = I\B of unblocked levels, with wavefunction 
ip(jl, ■ ■ ■ ,jn) (2j = Hjei\B denotes a sum over all unblocked levels in I). The cor- 
responding state l^n) is an eigenstate of the pair number operator and a Hamiltonian 
Hjj involving only pair operators: 

u 

5^6]6 i |* n )=n|* Tl ), Hu I *n)=£.|*n) , (13) 

3 

u 

A u = ^[2(e J - fi )6 ii - \d]b\b } . (14) 

ij 

Each eigenstate \^f n ,B) may be visualized as a coherent superposition of eigenstates 
of Hq that all lie in the same fixed-iV Hilbert space, and in all of which each pair of 
unblocked (j 6 U), time-reversed levels \j±) is either doubly occupied or empty. This 
is illustrated in Figs. ||(a) and (b), which schematically depict the exact ground states 
for even and odd N, respectively The odd ground state has a single blocked level, 
at the Fermi energy, containing an unpaired electron. The latter somewhat weakens 
pairing correlations relative to the even ground state and hence leads to parity effects, 
which will be extensively discussed in later sections. 

A useful measure for the amount of energy gained by \a) via its correlations is its 
"condensation energy" relative to the uncorrelated state |a)o, 

E c a ond =£ a - Q (a\H\a)o, where \a) = JJ c| ct JC/> , (15) 

and \U)o is the "Fermi ground state" in U , for which the n pairs occupy the n lowest- 
lying levels in U. 
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Note that Hjj is /i-independent, since the total Zeeman energy of any pair of 
electrons is zero. Hence the full /i-dependence of the eigenenergies resides in the 
rather trivial contribution 83(h) of the blocked levels, which is a very important and 
useful simplification. 

Diagonalizing Hjj would be trivial if the 6's were true bosons. However, they are 
not, and in the subspace spanned by the set U of all non-singly-occupicd levels instead 
satisfy the "hard-core boson" relations, 

b? = 0, [b p b\,] = 5 jj .{l-2tf j b j ), [b%b\,] = 5 jr b], (16) 

which reflect the Pauli principle for the fermions from which the 6's are constructed. 
In particular, fej 2 = implies that only those terms in (|ll]) arc non-zero for which the 
indices ji, ■ ■ ■ j n are all distinct. 

The task of finding the eigenstates \^ n ) is thus highly non-trivial. Nevertheless, 
an exact solution does exist. Unbeknownst to most of the condensed-matter physics 
community, it was found and studied extensively by Richardson in the 1960's and will 
be presented in Sec. [h]. Throughout the present part I, however, we shall use more 
well-known approaches based on the variational wavefunctions introduced by BCS , 
and that had been used to study the discrete BCS model before Richardson's solution 
was revived towards the end of 1998. 



5 Canonical characterization of pairing correlations 

Since the discrete BCS model has the standard reduced BCS form, the most natural 
first step toward understanding its T = properties is to use BCS-like variational 
wavefunctions (or equivalently Bogoljubov's mean-field approach), and indeed this 
will be done in Sees. |^ to |[ However, the limitations of such an approach should 
be realized at the outset: the spectra measured by RBT are excitation spectra for a 
grain with a fixed electron number N, and hence should be calculated for a grain with 
definite electron number N (i.e. completely isolated from the rest of the world, e.g. 
by infinitely thick oxide barriers). In contrast, the variational wavefunctions of BCS 
[Eq. (0) below] do not have the fixed-iV form [Eq. (|TT1)1 which any true eigenstate 
should have, but instead are formulated in a grand- canonical (g.c.) framework (as is 
the Bogoljubov mean-field approach to which they are intimately related). 

When considering a truly isolated superconductor such as a perfectly insulated 
grain (another example would be a superconductor levitating in a magnetic field due 
to the Meissner effect), one therefore needs to address the following question, which 
will be the main theme of the present section: how is one to incorporate the fixed-N 
condition into BCS theory, and how important is it to do so? This issue is well un- 
derstood and was discussed at length in the early days of BCS theory (Rickayzen's 
book gives a beautiful discussion [f75f ), in particular in its application to pairing cor- 
relations in nuclei |76|, p. 439] (see also the general remarks in j77j). Nevertheless, 
for pedagogical reasons the arguments are well worth recapitulating in the present 
context. 

We shall first remind the reader that the use of a g.c. framework is only a matter 
of convenience, since the essence of the pairing correlations that lie at the heart of 
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BCS theory is by no means inherently g.c. and can easily be formulated in canonical 



language |17 . We then show how standard BCS theory fits into this scheme, point out 
that the differences between results obtained using g.c. and canonical wavefunctions are 
negligible for d <C A, and conclude that for the purposes of gaining a phenomenological 
understanding of the experimental data, standard grand-canonical BCS theory should 
be sufficient. Nevertheless, the fundamental question of how to improve on this theory, 
in order to achieve a truly canonical description and to properly treat fluctuation 
effects, which become important for d > A [^8[ |7{| ^3, is interesting and important 
in its own right and will be addressed at length in Sec. 9.4. 

For simplicity, throughout the present subsection ^ we shall consider only the even 
ground state in the thermodynamic limit (in which even-odd differences are negligible), 
so that U = I and blocking effects need not be worried about. 



5. 1 The grand-canonical BCS wavefunction 

Conventional BCS theory describes the pairing correlations induced by an attractive 
pairing interaction such as -ff ro d of Eq. (|^) within a g.c. ensemble, formulated on a Fock 
space of states in which the total particle number N is not fixed. This is illustrated 
by BCS's famous variational ground state Ansatz 

|BCS) = Yiiuj + e^Vjb]) |Vac) , with u) + v) = 1, (17) 

3 

where the variational parameters Uj and Vj are real and <j)j is a phase (which, it turns 
out, must be ^-independent, for reasons discussed below). |BCS) is not an eigenstate of 
N and its particle number is fixed only on the average by the condition (A^bcs = N, 
which determines the g.c. chemical potential /i. Likewise, the commonly used g.c. 
definition 

A gc = XdJ2( b i)BCS = Xd^ujVjj+i (18) 

3 3 

for the superconducting pairing parameter only makes sense in a g.c. ensemble, since 
(bj) would trivially give zero when evaluated in a canonical ensemble, formulated on 
a strictly fixed- N Hilbert space of states. (We shall use the term "pairing parameter" 
instead of "order parameter" , since the latter carries the connotation of a phase tran- 
sition, which would require the thermodynamic limit N — > oo, which is not applicable 
for ultrasmall grains). 



5.2 Canonically meaningful definition for the pairing parameter 

A theory of strictly fixed- N superconductivity must evidently entail modifications of 
conventional BCS theory. However, these are only of technical, not of conceptual 
nature, since the essence of the pairing correlations discovered by BCS can easily 
be formulated in a canonically meaningful way, including a definition for the pairing 
parameter. We shall now attempt to explain, in intuitive, non-technical terms, how 
this may be done (our discussion is indebted to that of Rickayzen |7f|). Readers with a 
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preference for rigor may consult Sees. (KM to 12 for a corroboration, using Richardson's 



exact solution, of the arguments presented below. 

Let |G) be the exact even ground state of the system, depicted in Fig. ||(a). As 



explained in Sec. 4.4, it is a coherent superposition of eigenstates of Hq that all have 
the same N and in all of which each pair of time-reversed levels |j±) is either doubly 
occupied or empty. Due to this coherent superposition, |G) entails strong pairing 
correlations, whose essential properties may be understood by investigating how they 
modify the correlators 

dj = (6t6.) , rj C n (&t 6j ) , u) ee (6 3 .ftt) , (19) 
relative to the form these take on for the Fermi ground state |Fjv): 

(C ii )F = *«(^)F, (^) F = fl(-6i), (u 2 j ) F = 9(e j ). (20) 
Cij(— Cl-j) is the matrix element for the interaction to be able to scatter a pair of 
electrons from level j to i, and u| and are the probabilities to find level j doubly 
occupied or empty, respectively. The pairing correlations in |G) must be such that 
iJrod lowers the ground state energy below that of the uncorrelated Fermi sea |Fjv) by 
an amount that is extensive (oc N cx d^ 1 ) in the thermodynamic limit. Clearly, this 
requires that (H rc d)G ~~ (Hred)F is negative and extensive, i.e. that 

Xd^2 [ C U - - Xd J2J2 2Re ( C y) K N ( and positive). (21) 

ij i j<i 

In the second expression we neglected the diagonal terms, since their number is so small 
(only oc N) that Xd^jl^j ~ (^)f] is at best of order unity in the thermodynamic 
limit. For Eq. ( |2l| ) to hold, |G) must have two properties: 

(i) the number of Cy 's that differ significantly from zero (i.e. are of order unity) 



should scale like N 2 , i.e. one power of N per index 75 p. 167] 



(ii) most or all of the dj for i < j should have the same phase, since a sum over 
random phases would average out to zero. 

Since a suitable pairing parameter should vanish in the thermodynamic limit unless 
both these conditions hold, the definition 

A L„ = (Ad) 2 J2(Cij - (4 + c j+ ){4_ Cj _)) (22) 



(or its square root) suggests itself, where the subscript emphasizes that (22) is mean- 
ingful in a canonical ensemble too, and we subtracted^ the "normal-state contribution 

g This subtraction was suggested to us by Moshe Schechter, who pointed out that then Eq. ( ^ ) has 
a natural generalization to position space: it is the spatial average, A 2 an = (Ad) 2 /dfi dr*2 T(fi , T2 ) , 
of the two-point function 

T{r x ,r 2 ) = {ip\_(ri)i> 1 l(r l )tp_(r2)ip + (r2)) - (ifl\_(fi)ip + (rb)) W>1 (n)i/>_ (f 2 )) 

(with iperif) = Vol -1 / 2 52 £ e lfc r cg CT ), which evidently measures the amplitude for the propagation of 
pairs as opposed to uncorrelated electrons. Other definitions for a canonically meaningful pairing pa- 
rameter have been suggested |3i],[hJ, such as AdJ^j UjVj or Ad 52j K^l^j) ~ ( c j+ c j+) ( c j — c j — ^ 1 ^ 2 > 
but these focus only on requirement (i) and fail to incorporate requirement (ii). A quantity very sim- 
ilar to Eq. |22| was recently proposed in Eq. (55) of Ref. |26| , namely (Ad) 2 52ij [Cy — (CVj)a=o]- 
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to Cij." If (i) and (ii) hold, A can will take on a finite value; its relation to a gap in the 
spectrum will become clear below. In the bulk limit, A can can be shown [see Sec. £1 
to reduce to the "bulk pairing parameter" A of Eq. (||) . 

5.3 Redistribution of occupation probability across £f 

Now, property (i) can be realized if all Cij in a finite (d-independent) range of Ej's 
and Ej'a around the Fermi energy sp differ significantly from zero; the width of this 
range will evidently determine the magnitude of A can (provided (ii) also holds) , which 
conversely can be viewed as a measure of this width. But a nonzero Cij evidently 
requires both b\bAG) ^ 0, implying (%)g 7^ and (u;)g 0, and also (G|fo|fo ■ ^ 0, 
implying (v,)g 7^ and (%)g 7^ 0. The product (%Uj)g must thus be different from 
zero [in contrast to (ujVj)p = 0] for all Ej within a finite range around £p (cf. Fig. [?]). 
This can be achieved by smearing out the sharp steps of the 0-functions of (%)f and 
(mj)f, so that (vj)o [or (%)g] is nonzero also for a finite range of Ej above [or below] £p. 
In other words, for |G) some occupation probability must be redistributed (relative to 
|Fjv)) from below to above £f, as illustrated in Fig. ||. This redistribution, which was 
called pair- mixing in [p^ , |l9f , frees up phase space for pair scattering and so achieves 
a gain in interaction energy (provided (ii) also holds) that more than compensates for 
the kinetic energy cost incurred thereby. 

Furthermore, note that properties (i) and (ii) imply, even without detailed calcu- 
lations, that the spectrum will be gapped. Consider, for example, a "blocking excita- 
tion" that disrupts pairing correlations by having \j+) definitely occupied and \j—) 
definitely empty; since pair-scattering involving level j is blocked, the energy cost is 

(Ej - M ) - [(Ej - M)2(6t6 i )) -XdJ2 (b\bj + &}&,.)] (23) 

<(#) 

= (Ej - M )(l - 2^ 2 ) + A d {Cij + Cji) , (24) 

in which the restriction on the sum reflects the blocking of scattering events involving 
level j. Since the first term of ( p4[ ) is positive definite (particle-hole symmetry ensures 
that (^ — Vj) < if Ej — (j. < 0) and the second of order A can , the excitation energy 
will be finite even for d — > 0, implying the existence of a gap of order A can . Similarly, 
"phase-breaking excitations" that violate the fixed-phase condition (ii) are gapped too: 
for example, if (CV,) cxci tcd = Oground for a given j and all j), the energy cost 
is -AdX)i(^j)[(Cy + C'j* l )cxcited - (Cy + CjO ground] , which is at least of order 2A ca „. 

We see, therefore, that the essence of pairing correlations can readily be formulated 
in a canonical framework: (i) a redistribution of occupation probability across £f 
occurs, such that each level j in a finite range around £f has a finite probability of 
both being doubly occupied or empty, and (ii) any two components of the ground state 
wavefunction that differ only by the exchange of a pair of electrons between two levels 
i and j have the same phase. 

Since pairing correlations with these properties are the microscopic property at 
the heart of all manifestations of "superconductivity" , it seems reasonable to call a 
sample "superconducting" as long as it exhibits pairing correlations with measurable 
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Fig. 7 A cartoon depiction of why "superconductivity breaks down" when the sample be- 
comes sufficiently small. Vertical lines are drawn at each single-particle energy ej, spaced 
with a mean level spacing d corresponding to (a) a "large" grain (d <C A); (b) a "small" grain 
(d ~ 0.25A); (c) an "ultrasmall" grain (d ~ A). The lines' height represents the function 
UjVj = j A 2 /(e 2 + A 2 ) of standard bulk BCS theory, to illustrate the energy range (of width 
A around ep) within which pairing correlations are strongest. Loosely speaking, the number 
of single-particle levels in this regime, namely A/d, corresponds to "the number of Cooper 
pairs" of the system. Evidently, this number becomes less than one when d > A as in (c), so 
that it then no longer makes sense to call the system "superconducting" [cf. Sec. pjj]. 



consequences. And by this criterion the gap observed in the even grains of RBT 
certainly qualifies. 



5.4 Gauge symmetry breaking 

Note that property (ii) will be preserved under the gauge transformation Cj a —>■ e'^Cjo-, 
i.e. Cij — > e _2l ^ i_ ^-'Cy, only if all 0'- are equal, say (f>'j = </>' . Property (ii), and 
likewise the pairing parameter A can , therefore (a) are not gauge invariant "locally" in 
j-space, but (b) are gauge invariant globally. These are obvious consequences of the 
facts that (a) a correlated fixed-N state consists of a phase- coherent superposition of 
many different components, and hence cannot be invariant under arbitrary changes of 
the phases of individual components; and that (b) all of these components contain the 
same number of electrons N and hence under a global gauge transformation all pick 
up the same phase factor e lN ^ . Obviously, global gauge symmetry can therefore never 
be broken in a canonical ensemble. In contrast, the breaking of global gauge symmetry 
by the g.c. pairing parameter A gc of Eq. fll8|), which transforms as A gc — > e 2l< ^ A gc , is 
an inevitable consequence or artefact of its g.c. definition |7f| p. 142]. 
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5.5 Making contact with standard BCS theory- 
One of the breakthrough achievements of BCS was, of course, to propose a simple 
variational ground state which has precisely the properties (i) and (ii) described above: 
when evaluating the correlators of Eq. (|19|) using |BCS) of Eq. (|17|), one finds 

(%)lcs = u r fe)ics = v j> (Cu)bcs = u i v i u j v j e~ l ^ i ~' l ' i) , (25) 

and also (A^ an )ecs = |A gc | 2 - The definite-phase requirement (ii) can thus be imple- 
mented by choosing all the phases <f>j to be the same, say <fij = <fi for all j, thereby 
breaking local gauge invariance (usually one simply takes = 0); and requirement (i) 
is fulfilled automatically when minimizing the expectation value (-ff)BCS w.r.t. Uj and 
Vj, since this does yield smeared-out step functions, namely jl], |58| 



v] = i [1 - ( Ej - n)/Ej] , Ei= ^(e.-Zi) 2 + |A gc p. (26) 

Here we neglected terms that vanish for d —* 0, and A gc is determined by the famous 
gap equation (for T = 0), 



£j<UD 



(27) 



The BCS wavefunction instructively illustrates some of the general properties discussed 
above. Firstly, the product Ujv'j, shown in Fig. ^, has a bell-shaped form with a well- 
developed peak around £f of width ~ |A gc |, illustrating that pairing correlations are 
strongest within a region of width | A gc | around the Fermi surface. Secondly, the energy 
of a blocking excitation [Eq. fl24|)] reduces to (sj — u)(l — 2Vj) + 2ujVj\A gc \ = Ej, which 

is just the well-known energy of the Bogoljubov quasiparticle state 7j + |BCS), where 

7jCT = UjCj a - crvje^c]^ . (28) 
Thirdly, an example of a phase-breaking excitation is 

7 ] +7 j_ |BCS) - (-tye" 1 * + Uj b]) J] ( Ui + Vj J+bl)\Vac) , (29) 

which has (Cy) exc ited = —UiViUjVj and energy 2Ej. 

It should be appreciated, however, that BCS chose a grand- canonical construction 
purely for calculational convenience (as is made clear on p. 1180 of their original paper 
|Q): the trick of using a factorized form of commuting products in (|l7|), at the cost of 
TV-indefiniteness, makes it brilliantly easy to determine the variational parameters Uj 
and Vj. In fact, BCS proposed themselves to use the projection of |BCS) to fixed N 
as the actual variational ground state, namely J75| 

,-2k 



/>Z7T 

|PBCS) = / d^e-^TT^+e^jfctjIVac) (30) 
Jo 

1 

(JV/2) 
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(PBCS for Projected BCS), which is of the general form of Eq. (|TT|) . In the bulk 
limit (d/A <C 1), however, it is completely adequate to use |BCS): firstly, the relative 
error which its factorized form causes, by taking the occupation amplitude of level j 
to be independent of that of level i, scales like 1/N pp. 150,163]; and secondly, the 
fluctuations in its particle number, (A./V 2 )bcs = (A 2 )bcs — A 2 — J2j(^ u j v j) 2 ; are 
equal to irA/d in the bulk limit, in which the relative fluctuations (AA 2 )bcs/A 2 oc 
dA/eF 2 therefore vanish. Thus, bulk results obtained from |PBCS) or |BCS) are 
essentially identical. In fact, Braun jH], |42) checked by explicit calculation that the 
functions (vJ)g, (w 2 )pbcs and v 2 are practically indistinguishable even for d/A as 
large as 0.5 [see Sec. [l^]. Significant differences do develop between them once d/A 
increases past 0.5, however, as will be discussed in Sec. [l2|. 

To end this section, note that Fig. offers a very simple intuitive picture for why 
pairing correlations weaken with increasing level spacing until, in Anderson's words 
|yj , "superconductivity is no longer is possible" when d > A: an increase in level 
spacing implies a decrease in the number of levels within A of ef for which w 2 v 2 
differs significantly from zero, i.e. a decrease in the number of pairs with significant 
pairing correlations. This number, namely A/d, can roughly speaking be viewed as 
the "number of Cooper pairs" of the system, and when it becomes less than one, as in 
Fig. 0(c), it no longer makes sense to call the system "superconducting" . However, this 
should not be taken to imply that pairing correlations cease altogether in this regime; 
remnants of them do persist, in the form of fluctuations, up to arbitrarily large d/A, 
as will be discussed in detail in Sec. [12. 



6 Generalized variational BCS approach 

In the next several sections we review the generalized variational BCS approach used by 
Braun et al. j3l], [ll| [l2| to describe the paramagnetic breakdown of superconductivity 
in nm-scale grains in a magnetic field. This theory produces theoretical excitation 
spectra that are in good qualitative agreement with the measurements of BRT shown 
in Fig. ^| and thereby yields the most direct confirmation available of the relevance to 
experiment of the discrete BCS model. Moreover, it sheds considerable light on how 
"superconductivity breaks down" (more precisely, how pairing correlations weaken) 
with increasing d and h: As mentioned in the previous paragraph, in grains with 
d ~ A (bulk gap), near the lower size limit |^l[ of observable superconductivity, the 
number of free-electron states with strong pairing correlations (those within A of e-p) 
is of order one. Thus, even in grains in which a spectral gap can still be observed, 
pairing correlations are expected to become so weak that they might be destroyed by 
the presence of a single unpaired electron [fl7| . This can be probed directly by turning 
on a magnetic field, since its Zeeman energy favors paramagnetic states with nonzero 
total spin. 

The theory reviewed below exploits analogies to thin films in a parallel magnetic 
field |6^, [37), but explicitly takes account of the discreteness of the grain's spectrum. 
Since in RBT's experiments the temperature T = 50mK is much smaller than all other 
energy scales (d, A), we shall neglect finite-temperature effects and set T = 0. In 



Sees. 6.1 the eigenenergies £ a of the grain's lowest-lying eigenstates |a) are calculated 
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\s=i,k=2) 



state ||); (c) the spin-| ground state 
levels are drawn for h = 



Fig. 8 Cartoon of four typical variational states, labeled using the notation of Eq. (|33| ) for 
(a-c) and Eq. ( ^ ) for (d). They represent (a) the even ground state jO); (b) the odd ground 

I); (d) a spin- 1 excited state ||, 2). The single-particle 
with the chemical potential half-way between levels and 1 for 
even systems (a), but exactly on level for odd ones (b,c,d). The ellipses joining states on 
the same level are meant to represent a "Cooper pair" , and signify its being empty or doubly 
occupied with amplitude (uj + i)ji>j); solid (dashed) ellipses are used for levels that would be 
completely filled (empty) in the absence of pairing correlations. 



approximately using a generalized g.c. variational BCS approach that goes beyond 
standard mean- field theory by using a different pairing parameter A Q for each \a). 
The £ a are then used to discuss various observable quantities, such as ^.-dependent 
excitation spectra (Sec. ||) and various parity effects (Sec. ||). 

The reasons for deciding to calculate the excitation spectra, despite their fixed 
N nature, within a grand- canonical framework are as follows: Firstly, its simplicity. 
Secondly and perhaps most importantly, the exact eigenenergies have the general 
form E a = £ n + £b(/i) [Eq. (|r2|)], in which all h-dependence resides in the exactly 
known contribution £b(/i) from the blocked levels. The choice of approximation scheme 
therefore only affects £ n , which determines the h — properties of the spectrum, 
such as the size of the zero-field spectral gap, etc., but not the qualitative features 
of the /i-dependence. In particular, this means that all of the analysis below could 
easily be "made exact" by simply replacing the g.c. approximations for £ n by the 
exact values from Richardson's solution. However, this is expected to cause only 
slight quantitative differences, since, thirdly, canonical calculations (mentioned after 
Eq. ( pi] ) and discussed in Sec. 9.4) yield very similar results to g.c. ones as long as 
d/A< 0.5, which, by inspection of Fig. |[ does seem to be the case for the grain in 
question (the analysis of Sec. || yields d/A ~ 0.67). 



6.1 The generalized variational Ansatz 

The Zeeman term in the Hamiltonian of Eq. (^]) favors states with a nonzero total 
z-component of the total spin, s = ^J2ja ac ]a c ja (henceforth simply called "spin"). 
Increasing h will thus eventually lead to a series of ground state changes to states with 
successively larger spins. In general, therefore, we are interested in pair-correlated 
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states with nonzero spin, and in particular in their eigenenergies. Following Braun 
et al. jl8|, [l9], ^|), we now show how this can be calculated variationally, using the 
following general BCS Ansatz for a state \s, B) with N = 2n + 2s electrons and a 
definite total spin s > (first introduced by Soloviev for application in nuclei ]74|): 

i s ' b ) - n c i rM' s) + v J s,B) fy i yac )- ( 32 ) 

If the spin is nonzero, it is built up by placing 2s unpaired spin-up electrons in a set B 
of b = 2s single-particle levels [cf. Eq. while the remaining single-particle levels 

have BCS-like amplitudes to be either empty (u^' 3 ^) or doubly occupied by a pair 
(v^ B) ), with {uf B) ) 2 + (v^ B) ) 2 = 1. The subscript U over products (and over sums 
below) indicates exclusion of the singly occupied levels in B, for which vS s,B \ v^ s,B ^ 
are not defined. The product Y[j thus constitutes a g.c. approximation to the state 
|* n ) of Eq. (|p). 

More specifically, in a given spin-s sector of Hilbert space the following two types 
of specializations of Eq. (132) were studied in detail (p = 2s mod 2): 



s-l+p/2 U 

II illK + ^W' ( 33 ) 

i=-s+p/2 j 

\s,k) = c[ s _ 1+p/2+fc)+ c (s _i +p/2)+ |s) . (34) 

|s) is the spin-s state with the lowest energy, i.e. the "variational spin-s ground state", 
obtained by placing the 2s unpaired electrons as close as possible to £f [Fig. ||(b,c)], in 
order to minimize the kinetic energy cost of having more spin ups than downs. |s, k) 
is a particular type of excited spin-s state, obtained from |s) by moving one electron 
from its topmost occupied level (s — 1 + p/2) upwards by k units of d into a higher 
level (s — 1 +p/2 + k). These constructions are illustrated in Fig. || of which (a) and 
(b) represent the variational ground states of a grain with an even or odd number of 
electrons, respectively. 

The orthogonality of the wavefunctions, (s, B\s' , B') — 5 ss >8bb', implies that the 
variational parameters v^' 3 ^ and u^ s,B ^ must be found anew for each (s, B) (hence 
the superscript), by minimizing the variational "eigenenergies" 

£^ s (h,d) = (s,B\H\s,B) (35) 
= -2sh + J> 4 [2( 6j - ri(v^ B) ) 2 + Ad(# S) ) 4 ] (36) 



a 4e 



ieB j 
u „ 

(s,B) (s,B) N 



which we use as approximations to the exact eigenenergies £^ ct (h, d). The vj term 
is not extensive and hence neglected in the bulk case where only effects proportional 
to the system volume are of interest. Here it is retained, since in ultrasmall systems it 
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is non-negligible (but not dominant either) Il8l [Hi . Solving the variational conditions 

d£ BCS 

— fr%Y = in standard BCS fashion yields 

(„(••*))* = |(1 - + Al B }^) , & ee at /i - \d(v^f , (37) 

where the "pairing parameter" A sB is defined by the relation 

A S , B ee Adf>^ s ^, or i = dg 1 2 , (38) 

which in the limit d/A — ► reduces to the standard bulk T = gap equation. Note 
that it is /i-independent, because it involves only unblocked levels j £ U , which are 
populated by pairs with zero total Zeeman energy. Note also that in Eq. d37j ) the 
Xd(v < j S ' B ' 1 ) 2 shift in £j, usually neglected because it simply renormalizes the bare ener- 
gies, is retained, since for large d it somewhat increases the effective level spacing near 
£f (and its neglect turns out to produce a significant upward shift in the £f B s (h, d)'s, 
which one is trying to minimize). 

The chemical potential /i is fixed by requiring that 

u 

2n + 2s = (s,B\N\s,B) = 2s + 2^(^. S ' B) ) 2 . (39) 

3 

In contrast to conventional BCS theory, the pairing parameter A s ^b can in general not 
be interpreted as an energy gap and is not an observable. It should be viewed simply 
as a mathematical auxiliary quantity which was introduced to conveniently solve the 
variational conditions. However, by parameterizing v^'^ and A s- b does serve 

as a measure of the pairing correlations present in \s, B): for vanishing A S: b the latter 
reduces to an uncorrelated paramagnetic state \s,B)q with spin s and energy £® B , 
namely 

u 

\s,B) = II 4 + n &]|0> , with £° <B = (s,B\H\s,B) , (40) 
ieB j<a 

and the condensation energy E^° B d = £f B s — £® B of \s, B) reduces to zero. 



6.2 General numerical solution - illustration of the blocking effect 



The simultaneous solution of Eqs. (|37|), (|38|) and (|39j) is a straightforward numerical 
exercise which Braun and von Delft performed 18, [19], for the sake of "numerical con- 
sistency", without further approximations. (Analytical solutions can be found only in 
the limits d<C A and d 3> A Si b, see App. A of |l!|.) The numerical results are sum- 
marized in Fig. ^, which shows the pairing parameters A s _b [Figs. ^(a,c)] and energies 
£® B S [Figs. ^|(b,d), solid lines] of some selected variational states \s,B), as well as 
the energies £® B of the corresponding uncorrelated states \s, B)q [Figs. ^(b,d), dashed 



Jan von Delft, Superconductivity in ultrasmall metallic grains 



27 




0.0 0.5 1.0 0.0 1.0 2.0 

Level spacing d/A Level spacing d/A 

Fig. 9 Properties of (a,b) spin-s ground states |s) = \s, 0) [Eq. ([sif)] and (c,d) spin-i 
excited states ||,fe) (for fc = 0, ...,4) [Eq. (|3J)], as functions of d/A (i.e. decreasing grain 
size), calculated for A = 0.194. (a) The pairing parameters A s (d)/A, which vanish at critical 
level spacings d? cs of 2.36, 0.77, 0.44, 0.31,. . . for s = 0, 1/2, 1, 3/2, . . ., respectively, (c) The 
pairing parameters A 1 ^ 2 ,k, together with their k — > oo limit, which equals Aq (dashed line). 
(b,d) show relative energy densities (since normalized by d/A 2 oc Vol -1 ) at h = for both 
correlated (solid) and uncorrelated (dashed) states, the latter obtained by setting A S: b = 
in the former, (b) shows (£f cs -£ ° /2 )d/A 2 (solid) and (£° ~£° /2 )d/A 2 (dashed), the energy 
differences of \s) and |s)o relative to the uncorrelated spin-p/2 Fermi sea \p/2)o- (d) shows 
{£ff 2 s k -£<l /2 )d/A 2 (solid) and (£° /%k - £° /2 )d/A 2 (dashed), the energy differences of 
and ||,fc)o relative to the uncorrelated spin-^ ground state ||,0)o = | )o - Solid and dashed 
lines meet at the critical level spacing dfjl s at which A Sj ft becomes and the condensation 
energy E c a J d = £f^ s - £° e>k vanishes. 



lines] ; both £fg and £ ° B are plotted relative to the energy £^ 2 of the uncorrelated 
spin-p Fermi sea |p/2) . The results have a number of salient features: 

(i) In the bulk limit d/A — > 0, all of the pairing parameters A Si b reduce to A, as 
expected, and the energy differences 

£?3 S ~ £ p/2 between the correlated states \s,B) 
and the uncorrelated Fermi sea \p/2) reduce to —\A 2 /d = —^N{e-p)^ 2 , which is the 
standard bulk result for the condensation energy. 

(ii) Each A Sj # in Figs. ^(a,c) decreases with increasing d. This reflects the fact that 
with increasing d, the number of pair-correlated states within A of e-p decreases [cf. 
Fig. |^ and the last paragraph of Sec. 5.5 1, so that the amount of pairing correlations, 
for which A s ^b is a measure, decreases too. 

(iii) Each A Sy B vanishes abruptly at a critical level spacing df^ (whose precise 
numerical value depends sensitively on model assumptions such as the value of A and 
the use of uniformly-spaced levels |2(|). For d > <i^ s no pairing correlations exist 
at this level of approximation, so that that the condensation energy E^°g d (difference 
between solid and dashed lines) vanishes and the solid and dashed lines in Figs. §(b,d) 
meet. 

(iv) In Figs. ||(a), the pairing parameters A s for the spin-s ground states decrease 
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rapidly with increasing s at fixed d (and df cs < df, cs if s > s'). [This is a generaliza- 
tion of a parity effect discussed by von Delft et al. |l7j] , who studied only ground state 
pairing correlations and found that these are weaker in odd (s = 1/2) grains than in 
even (s = 0) grains, A1/2 < Ao, cf. Sec. @.] This tendency is a direct consequence of 
the blocking effect described in Sec. and is independent of model details: larger s 
means more unpaired electrons, more terms missing from the sum ■ , less correlated 
pairs and hence smaller A Sv b. 

(v) As d increases the blocking effect described in (iv) becomes stronger, i.e. the 
difference between the various A s for different s becomes more pronounced, since then 
the relative weight of each term missing in the sum yj^ increases. The blocking effect 

is most dramatic in the regime d/A £ [0.77, 2.36] in which Ao ^ but A^o = 0. 
This is a regime of "minimal superconductivity" , in the sense that all pairing 

correlations that still exist in the even variational ground state |0) (since Ao ^ 0) are 
completely destroyed by the addition of a single electron or the flipping of a single 
spin (since A s ^ = 0). 

(vi) Considering the spin-i excited states ||, k) of Figs. ^(b,d), one finds that the 
larger k, the longer the pairing correlations survive with increasing d: the critical 
spacings d®^ k increase with k, approaching the value d^ cs of the spin-0 case as 
k — ► oo; correspondingly, the larger k, the larger the c£-value at which the condensation 
energies Ey^ d k [differences between solid and dashed lines in Fig. ^(d)] vanish. The 
intuitive reason why the amount of pairing correlations in an excited \s, k) increases 
with k is of course quite simple: the further the unpaired electron sits from the Fermi 
surface where pairing correlations are strongest, the less it disrupts the latter (since 
uuVk becomes very small for large k, see Fig. 0). In fact, the state \\,k — > oo) will 
have just about the same amount of pairing correlations as the even ground state |0) 

(Ai/2,fc->oo ^ A ). 

(vii) Similar effects hold for excited states in other spin sectors (not shown): The 
higher the excitation, the larger the pairing parameter A s ^b- However, the concomit- 
tant gain in correlation energy is always less than the kinetic-energy cost of having an 
unpaired electron far from ep. 

(viii) The strong dependence of A Sj s on s and d for d > A illustrates why in this 
regime a conventional mean-field treatment is no longer sufficient: the system cannot be 
characterized by a single pairing parameter, since the amount of pairing correlations 
vary from state to state, each of which is characterized by its own pairing parameter. 
Instead, the present variational approach is, roughly speaking, equivalent to a doing 
a separate mean-field calculation for each new choice U of unblocked levels within 
the Fock space spanned by them (i.e. replacing bj — > {bj — (bj)} + bj and neglecting 
terms quadratic in the fluctuations {bj — (bj)})- Indeed, the behavior of A s ^(d) near 

dfjp has the standard mean-field form ^/l — d/df^p, as can be shown analytically 

IB A PP . A]. 

To summarize: pairing correlations decrease with increasing d and s and decreasing 
k. These features survive also in more accurate canonical calculations. This is not the 
case, however, for the abrupt vanishing of A s b at df^ , which signals the breakdown of 

the g.c. approach once d becomes of order A: canonical methods show that, regardless 
how large d becomes, some remnants of pairing correlations survive and the pairing 
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parameters (A s .s) C an do not vanish [Sec. |9.4| ], in accordance with the rule of thumb 
that "in a finite system no abrupt phase transition can occur between a zero and 
nonzero order parameter." 



7 Softening of the iJ-induced transition to a paramagnetic state 

Since states with nonzero spin are favored by the Zccman energy but have smaller 
correlation energy due to the blocking effect, a competition arises between Zeeman 
energy and correlation energy. The manifestations of the blocking effect can thus be 
probed by turning on a magnetic field; if it becomes large enough to enforce a large 
spin, excessive blocking will destroy all pairing correlations. 

The situation is analogous to ultra-thin films in a parallel magnetic field 1 66 , ^| , 



where orbital diamagnetism is negligible for geometrical reasons and superconductivity 
is destroyed at sufficiently large h by Pauli paramagnetism. This occurs via a first or- 
der transition to a paramagnetic state, as predicted by Clogston and Chandrasekhar 
(CC) |l) by the following argument (for bulk systems): A pure Pauli paramagnet 
chooses its spin s such that the sum of the kinetic and Zeeman energies, s 2 /Af(ep)— 2hs, 
is minimized, and hence has spin s = /iA/"(£f) and ground state energy — /i 2 A/"(£f)- 
When this energy drops below the bulk correlation energy — ^A 2 Af(eF) of the super- 
conducting ground state, which happens at the critical field hcc = A/^/2, a transition 
will occur from the superconducting to the paramagnetic ground state. The transition 
is first-order, since the change in spin, from to sec = ^ccA/"(£f) = A/(dV2), is 
macroscopically large (M(ef) — 1/d oc Vol). 

This transition has been directly observed by Meservey and Tedrow |]6| ^37) in 
ultra-thin (5nm) superconducting Al films (A = 0.38meV), whose density of states 
[Fig. [To|(a)l they measured via the tunnel conductance through an oxide layer between 
a normal metal and the film. They found that in a magnetic field the BCS quasiparticlc 
peak splits up into two subpeaks, separated in energy by 2^j, B H [Fig. |l^(b)], which 
simply reflects the Zeeman splitting of quasiparticles states^ with spin up or down 
(and g = 2). Remarkably, the tunneling threshold abruptly dropped to zero at a 
field of 4.7 T [Fig. [tO](b)], which they associated with the field Hqc at which the 
phase transition from the superconducting to the paramagnetic ground state occurs. 
Indeed, Fig. |l^(b) demonstrates clearly that the transition to the normal state is first 
order: the mean of the spin-up and spin-down spectral gaps, i.e. the pairing parameter 
A, is constant until the critical field Hqc is reached, at which it abruptly drops to zero. 

For the case of isolated ultrasmall grains, the above picture of the transition needs 
to be rethought in two respects due to the discreteness of the electronic spectrum: 
Firstly, the spin must be treated as a discrete (instead of continuous) variable, whose 
changes with increasing h can only take on (parity-conserving) integer values. Sec- 
ondly, one needs to consider more carefully the possibility of /i-induced transitions 
to nonzero spin states that are still pair- correlated (instead of being purely paramag- 
netic), such as the variational states \s,B) discussed above. (In the bulk case, it is 



h Recall that the BCS quasiparticles "/j a = u j c \ a ~ av j c j _ CT have well-defined spins 
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Fig. 10 Thin films in a magnetic field (Figs. 11 and 12 of |66|). (a) Tunneling conductance 
from a normal metal through a tunnel barrier into a thin superconducting Al film, as function 
of voltage, for several magnetic fields labeled in increasing order "a" to "f". The conductance 
reflects the BCS quasiparticle density of states, whose single peak (for a given sign of V) 
for H — splits into two separate peaks for H 7^ 0, corresponding to the Zeeman energy 
difference between quasiparticles with spin up and down, (b) Voltage corresponding to the 
maxima of spin-up and spin-down density of states as a function of magnetic field. At the 
critical field Hcc superconductivity is destroyed and the tunneling threshold drops abruptly 
to zero. 

obvious that such states play no role: the lowest pair-correlated state with nonzero spin 
obtainable from the ground state by spin flips is a two-quasiparticle state, costing en- 
ergy 2A — 2h; when h is increased from 0, the paramagnetic transition at hcc = A/v2 
thus occurs before a transition to this state, which would require h — A, can occur.) 

Quite generally, the effect of increasing h from can be analyzed as follows: At 
given d and h, the grain's ground state is the lowest-energy state among all possible 
spin-s ground states \s) having the correct parity p = 2s mod 2. Since £ s (h,d) = 
£ s (0,d) — 2/is, level crossings occur with increasing h, with £ a i dropping below £ s at 
the level crossing field 

t>>i(d) = M0g_|M (41) 

Therefore, as h is slowly turned on from zero with initial ground state \sq = p/2) , a cas- 
cade of successive ground-state changes (GSCs) to new ground states |si), \s2), ■ ■ ■ will 
occur at the fields h S0 S1 , h sl S2 , . . . Let us denote this cascade by (sq, Si); (si, S2); ■ ■ -i 
for each of its GSCs the corresponding level-crossing fields h SiS < (d) is shown in Fig. [I]] . 
Generalizing CC's critical field to nonzero d, let us denote the (parity-dependent) field 
at which the first transition (so,s±) occurs by hcc{d,p) = h SQ , Sl (d), which simply is 
the lower envelope of the level-crossing fields h So , Sl in Fig. [ll] (shown as bold solid and 
dashed lines for so — and so = 5, respectively). In the limit d — > it is numeri- 
cally found to reduce to the Clogston-Chandrasekhar value, i.e. hcc(0,p) = A/y/2, as 



In general, the order in which the GSCs occur with increasing h within a cascade 
(i.e. the order of h SlS i lines encountered when moving vertically upward in Fig. [ll]) 



expected. 
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Fig. 11 d-dependence of the level-crossing fields h SiS ,(d)/A [Eq. @] at which £f, cs drops 
below £f cs with increasing h. Only those level crossing fields are shown that belong to 
the cascade of (fixed-N) ground state changes (GSCs) (so,si); (si,S2); ...that occur as 
h increases from at given d. Solid (dashed) lines are used for even (odd) grains with 
integer (half-integer) spins, and some are labeled by the associated GSC (s,s r ). (In contrast, 
in Fig. [12] the N-changing tunneling transitions in are labeled by |sj) — ► |s/).) The size 
\A£ B1 j-i — A£ SOj /| of the first jump (occuring at the level-crossing field ftcc(p, d) = h so , ai ) in 
the lowest line of the tunneling spectra of Fig. |l^ is shown by the lowest two (jagged) curves 
(solid for e — ■> o and dashed for o —* e tunneling spectra), which both approach the CC value 
1 - 1/V2 as d -> 0. 

depends sensitively on d and an infinite number of distinct regimes (cascades) I, II, III, 
. . . can be distinguished: Starting at large d we find the typical normal-grain behavior 
(0, 1); (1, 2); (2, 3); . . . for even grains and (|, §); (§, §);••■ for odd grains, with ho,i < 
(or >) hi a in regimes I (or II). In regimes III and IV of somewhat smaller d, the 
order of GSCs is (0, 2); (2, 3); . . . and (3, §);(§,§);.. ., etc, i.e. the spin si attained 
after the first GSC (so,si) has increased to 2 in the even case. This illustrates a 
general trend: the spin si(d) after the first transition increases with decreasing d and 
becomes macroscopically large in the d — > limit, where s% = hcc/d = A/(dy/2), as 
explained in recounting CC's argument above. 

Furthermore, it turns out that A Sl (d) — and therefore £ = for all d, 
implying that after the first GSC the new ground state \s±) is always (not only in 
CC's bulk limit) a purely paramagnetic state, i.e. without any pairing correlations 
in the g.c. framework (canonical calculations would yield some weak remnant pairing 
correlations in the form of fluctuations). In this regard, CC's picture of the transition 
remains valid throughout as d is increased: at hcc(d,p), a transition occurs from 
the superconducting ground state to a paramagnetic, uncorrelated state |si)oj the 
transition being first-order in the sense that A Sl (d) — 0; however, the first-order 
transition is "softened" with increasing d, in the sense that the size of the spin change, 
si — sq, decreases from being macroscopically large in the bulk to being equal 1 at d^$> A 
(regimes I and II). 

To conclude this section, we mention that the above analysis of the paramagnetic 
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breakdown of superconductivity has recently been generalized to finite temperatures 
p8| , using the so-called static path approximation [explained in Sec. 14.3 1 to treat 
fluctuation effects properly. 



8 Excitation spectrum in a magnetic field 

In this section we compare the theoretical tunneling spectra for a grain coupled to 
leads, calculated as functions of h and d (l^, [lj| and compare these to RBT's 
measurements of Fig. 0. 

The form of the tunneling spectrum depends in a distinct way on the specific choice 
of level spacing d and on the electron number parity p of the final states |/) of the 
bottleneck tunneling processes \i) — > |/) (or \a') — ► \a) in the notation of Sec. |^). 
However, for the uniformly spaced £j-levels used here, particle-hole symmetry ensures 
that there is no difference between electron addition or removal spectra |ijvq=i) ~~ ¥ \In)- 
To calculate the spectrum for given d and p, Braun et al. []l8| , [l9| , [i"2|| proceeded 
as follows: they first analyzed at each magnetic field h which tunneling processes 
\i) — > |/) are possible, then calculated the corresponding tunneling energy thresholds 
A£ if {h) ee S f (h) - £i{h) [cf. Eq. (|)] and plotted A£ if {h) - A£ min (0) as functions of 
h for various combinations of \i) and |/), each of which gives a line in the spectrum. 
Since the selection rule Sf — Si — ±1/2 holds, only slopes of ±1 can occur. The 
reason for subtracting A£ m j n (0), the h = threshold energy cost for the first (lowest- 
lying) transition, is that in experiment, this energy depends on V g and hence yields 
no significant information (see also Sec. 2.4.6 of Ref. Jl6|]). Neglecting nonequilibrium 
effects H |[ [l(| (which were minimized in the present experiment by tuning Vg), 
the initial state is always taken to be the ground state of a given spin-s sector. The 
appropriate Si(h, d) must be determined from Fig. [if]. 

Fig. [l^ shows four typical examples of such theoretical tunneling spectra, with 
some lines labeled by the corresponding \i) — ► |/) transitions. Whenever h passes 
through one of the level-crossing fields h SuS ., of Eq. (|4l|), the grain experiences a 
ground state change (s^, Si<), at which the set of allowed tunneling transitions changes 
from \si) —> {I s /)} to | s^/) — > {|s/')l- Therefore, at /i Si , s ., one set of lines in the 
tunneling spectrum ends and another begins, producing kinks or discontinuities. A 
kink occurs if one of the new final states coincides with one of the old ones, \f) = \f), 
meaning that it can be reached from both \si) and \si>) [i.e. Sf — Si — — (sf — Sj/)], 
in which case A£if(h) and A£ifi(h) have slopes of opposite sign. However, for most 
lines this is not the case, so that at h s . s i the line \si) — > |/) simply ends while new 
lines | Si') — > \ f) begin. This results in discontinuities (or "jumps") in the spectrum 
at /i SijS ' of size (A£i>f> — A£if)(h s . s '), unless by chance some other final state \f) 
happens to exist for which this difference equals zero. 

Since the order in which the GSCs (s^, s^) occur as functions of increasing h depend 
on d and p, as indicated by the distinct regimes I, II, III, ... in Fig. |ll|, one finds a 
distinct kind of tunneling spectrum for each regime, differing from the others in the 
positions of its jumps and kinks. In regime I, where the order of occurrence of GSCs 
with increasing h is (0, 1); (^, |); (1, 2); (|, |); . . ., there are no discontinuities in the 
evolution of the lowest line [see Fig. |l2|(a)]. For example, for the e — > o spectrum, the 
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Fig. 12 The theoretical odd-to-even and even-to-odd tunneling spectra (A£if — A£ m i n (0))/A 
predicted for an ultrasmall superconducting grain as a function of magnetic field h, to be 
compared with the experimental spectra of Fig. |B|, for two different level spacings: (a) d = 
0.67A and (b) d — 0.34A (corresponding to regimes I and III of Fig. [ll| respectively). Some 
lines are labeled by the corresponding S{ — > s\ tunneling transitions. Not all possible higher 
lines (corresponding to excited final states \s,j)) are shown. Vertical dashed lines indicate 
those level-crossing fields h s<s i [Eq. (ffl|)1 at which kinks or jumps occur, with /to.i < fti/2,3/2 < 
hi,2 < h a / 2 ,B/2 in (a) and fei/2,3/2 < ho, 2 < ^2,3 in (b). 



lowest |0) — > 1 1 / 2} line changes continuously to |1) — > |l/2) at ho,i, since |s/ — s-| = 1/2. 
However, in all other regimes the first change in ground state spin (at /io, Sl from to 
si) is > 1, implying a jump (though possibly small) in all e — ► o lines, as illustrated 
by Fig. ||(b). 

The jump's magnitude for the tunneling thresholds, i.e. the lowest e — > o and o — ► e 
lines, is shown as function of in the lower part of Fi g. pd|. It starts at d = from 
the CC value A(l - 1/V$) measured for thin Al films |66|, |67[, and with increasing a 
decreases to (non-monotonically, due to the discrete spectrum) . This decrease of the 
size of the jump in the tunneling threshold reflects the fact, discussed in Sec. |7[ that 
the change in spin at the first ground state change (so, si) decreases with increasing d 
(as si — so ~ hec/d), and signals the softening of the first-order superconducting-to- 
paramagnetic transition. 

The fact that the measured tunneling thresholds in Fig. [5] show no jumps at all, 
which might at first seem surprising when contrasted to the threshold jumps seen at 
hec m Fi g. |10| fo r thin films in a parallel field @ |7j, can therefore naturally be 
explained Q |19|] b y assuming the grain to lie in the "minimal superconductivity" 
regime I of Fig. |ll| (where the jump size predicted in Fig. [ll] is zero). Indeed, the 



34 



Ann. Phys. (Leipzig) 10 (2001) 3 



overall evolution (i.e. order and position of kinks, etc.) of the lowest lines of Fig. 
qualitatively agrees with those of a regime I tunneling spectrum, Fig. ^^(a). This 
important result rather convincingly establishes the phenomenological success of the 
discrete BCS model. It also allows one to deduce the following values for the level- 
crossing fields H s . s i (indicated by vertical dashed lines in Figs. || and |l2|): Ho,i = 4T, 
#1/2,3/2 = 4.25T, iJi i2 = 5.25T and #3/2,5/2 — 6.5T. As corresponding uncertainties 
we take AH SiyS > = 0.13T, which is half the # resolution of 0.25T used in experiment. 

By combining the above H Si y values with Fig. [ll], some of the grain's less- well- 
known parameters can be determined somewhat more precisely: 

(i) To estimate the grain's "bulk #cc'\ note that since #i/2,3/2/#o,i — 1-06, 
this grain lies just to the right of the boundary between regions II and I in 
Fig. [ll] where d/A ~ 0.63, at which we have /io,i/^cc — 0.95, so that #cc = 
#0,1/0.95 ~ 4.2 T. This is quite close to the value H C c — 4.7 T found experi- 
mentally |66|, ^ in thin films in a parallel field, confirming our expectation that 
these correspond to the "bulk limit" of ultrasmall grains as far as paramagnetism 
is concerned. 

(ii) The grain's corresponding bulk gap is A = v / 2A i B#cc — 0-34 meV, implying 
a coupling constant of A = 0.189 [by Eq. (§)]. A posteriori, these values can 
be regarded as being more appropriate for the present grain than the choices 
A = 0.38 meV and A = 0.194 made in Sec. |T^, though the differences are 
evidently not significant (12% for A and 3% for A). 

(iii) The mean level spacing implied by d/A ~ 0.63 is d ~ 0.21 meV. The crude 
volume-based value d ~ 0.45 meV cited in the caption of Fig. g thus seems 
to have been an overestimate. It would be useful if this determination of d 
could be checked via an independent accurate experimental determination of d 
directly from the spacing of lines in the excitation spectrum. Regrettably, this 
is not possible: the measured levels are shifted together by pairing interactions, 
implying that their spacing does not reflect the mean independent-election level 
spacing d. Nevertheless, note that the measured spacing of 0.05 meV between 
the lowest two states of the odd grain agrees quite well with the crude BCS 
estimate VA 2 + d 2 -A [cf. Eq. @], which gives 0.06 meV when evaluated for 
d = 0.21 meV and A = 0.34 meV. 

The higher lines plotted in Fig. [l2] correspond to transitions into spin-s / state of 
the form \sf, k) [ cf . Eq. (|34| ) and Fig. ||(d)] (for simplicity these were the only ones 



considered in p8l 19, |42j, though in general others are expected to occur too). The 
jumps in these lines, e.g. in Fig. |l^(a) at /ii,2, occur whenever the two final excited 
states \sf, kf) and \sf>, kf) before and after the GSC at h Si>s '. have different correlation 
energies. (Recall that the correlation energy of an excited state \sf,Bf) can be nonzero 
even if that of the corresponding ground state \sf) is zero, since the former's unpaired 
electrons are further away from Sp, so that A Sf g f > A Sf , see point (vi) of Sec. |^.) 
Experimentally, these jumps have not been observed. This may be because up-moving 
resonances lose amplitude and are difficult to follow jjj with increasing h, or because 
the widths of the excited resonances (~ 0.13A) limit energy resolution §, g, |T|. 
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For somewhat larger grains, the present theory predicts jumps even in the lowest 
line, as illustrated in Fig. |l^(b). It remains to be investigated, though, whether orbital 
effects, which rapidly increase with the grain size, would not smooth out such jumps. 

To conclude this section, we emphasize once again that more than qualitative 
agreement between theory and experiment can not be expected, since both the model 
and our variational treatment thereof are very crude: the model neglects, for instance, 
fluctuations in level spacing and in pair-coupling constants, and the g.c. wave functions 
become unreliable for d/ A > 0.5. Furthermore, we neglected nonequilibrium effects in 
the tunneling process and assumed equal tunneling matrix elements for all processes. 
In reality, though, random variations of tunneling matrix elements could suppress some 
tunneling processes which would otherwise be expected theoretically. 



9 Measurable consequences of the blocking effect: parity effects 

This section is devoted to various measurable manifestations of the blocking effect, in 
the form of parity effects, i.e. differences between a grain with an even or odd number 
of electrons. 

9.1 Bulk consequences of blocking 

The most obvious measurable manifestation of the blocking effect is the very existence 
of a spectral gap: "breaking a pair" and placing the two newly unpaired electrons 
in two singly-occupied levels costs a significant amount of correlation energy, because 
the unpaired electrons loose pairing energy themselves and also disrupt the pairing 
correlations of the other pairs. This, of course, is already present in standard bulk 
mean-field BCS theory via the energy cost of at least 2A involved in creating two 
quasiparticles, and is one of the hallmarks of superconductivity. 

In the context of ultrasmall grains, let us denote the pair-breaking energies for an 
even (odd) grain, i.e. the minimum energy cost per electron for breaking a pair by 
flipping a single spin at h = 0, by Q, e (fl )' 

fi e = 5(^1 ~ £o)h=0, ^0 = 5(^3/2 — £l/2)h=0 ■ (42) 

The even pair-breaking gap f2 e is of course strikingly visible in RBT's h — spectra as 
a large spectral gap for even grains [cf. Figs. || and||; the latter gives fi e = 0.26 meV]. 
Its presence is direct evidence for the existence of pairing correlations in the grain, 
which in that sense can still be called "superconducting" . 

In contrast, the odd pair-breaking gap Q a can not be obtained from h — spectra, 
since in an odd grain the lowest excitation does not involve breaking a pair, but 
simply exciting the unpaired electron, which does not require a correlation-induced 
gap to be overcome. To measure Q a , a finite field is needed: by Eq. (|4l|), fl e — ft, .i 
and Q a — /11/2.3/2J hence both spin-flip gaps are equal to level-crossing fields that 
can be deduced from /i / data, as explained in Sec. ||. For Fig. [| this yields 
fl e = 0.23±0.01 meV and fi = 0.24±0.01 meV [a result further discussed in Sec. [OJ- 
The reason that the fi e -value determined in this way is somewhat smaller than the 
above-mentioned 0.26 meV determined at h = is presumably that the experimental 
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spectral lines are not perfectly linear in h (having a small /^-contribution due to orbital 
diamagnetism, which should cause the spectroscopic gap to decrease faster with h than 
in our model). 

Another consequence of the blocking effect is that the condensation energies £p° 2 d = 
£ p /2 — £® j2 f° r an even and odd grain differ: the unpaired electron of an odd grain weak- 
ens its pairing correlations relative to an even grain, so that Ey^ d is less negative than 
iS , Q° nd . In the bulk limit their difference approaches £ , 1 ™ d — Eg° nd — > A, the energy of a 

single quasiparticle. For large mesoscopic islands (with d/A -C 1) this energy difference 
has indeed been directly observed: it causes a change from e- to 2e-periodicity in the 
gate- voltage dependence of Coulomb oscillations (|^. |2[ |83|, |4[ ^5|, |8(| . For ultrasmall 
grains, however, ground state energy differences arc currently not directly measurable, 
due to experimental difficulties! explained in detail in Sec. 2.4.6 of Ref. 

The parity effects discussed above survive in the bulk limit. Let us now turn 
to parity effects that result from even-odd differences in the d- dependence of various 
quantities. 



9.2 Parity-dependent pairing parameters 

As is evident from Fig. ^(a,b), not only the condensation energies E p / 2 are parity 
dependent; as soon as one leaves the bulk regime, the pairing parameters A p / 2 become 
parity-dependent too, with A > Ai/g. In the context of ultrasmall grains this was 



first emphasized by von Delft et al. 17 , but it had been anticipated before by Janko, 
Smith and Ambegaokar |7j| and Golubev and Zaikin j8?|], who had studied the first 
correction to the bulk limit, finding Ao — Ajy 2 — d/2 to leading order in d/A; and 
this result, in turn, had already been published by Soloviev in the nuclear physics 
literature as long ago as 1961 Q . 

The g.c. results of Fig. ^(a), in particular the fact that the critical level spacing 
d^9 2 at which A p / 2 vanishes is smaller for odd than even grains {dfuf < d^ cs ), 
suggest that "pairing correlations break down sooner in odd than even grains" p7[ . 
However, it should be remembered that the vanishing of A p / 2 signals the breakdown 
of the g.c. approach. A mo re a ccurate statement, that is born out by the canonical 
calculations reviewed in Se c. |9.4| , is that the inequality £i / 2 > £ o persists for arbitrarily 
large d (see Fig. |lj in Sec. |ll|), i.e. pairing correlations are always weaker for odd than 
even grains, although they never vanish altogether in either. 



9.3 Matveev-Larkin parity parameter 

To capture the difference between correlations in even and odd grains in terms of 
measurable quantities (which A p / 2 are not), Matveev and Larki n |24|| proposed the 
parameter (sometimes called "pairing energy" in nuclear physics Jiff) 

A^ L = E^ 1 - ±(£^+£^+ 2 ) (where N is even) , (43) 



'In brief: when V g is changed over a sufficiently large range to see Coulomb oscillations (~ e/Cg ~ 
1 V), sudden rigid shifts in the background off-set charge Qo are encountered at random values of 
V g , which spoil the 2e-periodicity which would have been expected otherwise. 
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Fig. 13 The parity parameters (a) Ap L of Matveev-Larkin [Eq. @] and (b) A p b for the 
pair-breaking energies [Eq. (jig)], as functions of d/A, calculated using the g.c. variational 
BCS approach of Sec. S.l (dashed lines), and Richardson's exact solution of Sec. hd (solid 
lines). In (a), we also show the perturbative result for the uncorrelated Fermi sea, (A P ) pC rt = 
^Xd (straight dash-dotted line); and the renormalized result (A P 1L ) rcn ~ d/[2 ln(ad/A)] of 
Eq. (^J), in its range of validity d/A ^> 1 (dotted line). The parameter a = 1.35 is chosen to 
ensure quantitative agreement with the exact result in the limit d/A 2> 1. For a summary of 
the results of various other canonical calculations of Ap L , see Fig. [l4|. The inset of (a) shows 
the Dyson equation used to calculated the renormalized coupling A in Eq. (ft5|). 



i.e. the difference between the ground state energy of an odd grain and the mean of the 
ground state energies of the even grains obtained by removing or adding one electron. 



Fig. |l3|(a) shows its behavior as function of d/ A. In the bulk limit we have £$ 



£ 



N+2 



and Ap ~ A, which is simply the energy cost for having an unpaired electron on the 
odd grain. With increasing d/A, this energy cost decreases since pairing correlations 
get weaker, hence Ap L initially decreases. It begins to increase again for d > A, since 
then pairing correlations are so weak that the behavior of is governed by the 



"self-energy" of the one extra pair in £q +2 relative to £q . For example, in the g.c 



variational BCS result for Ap L , 



l^p )bcs 
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namely 

■Ad/2, 



(44) 



it is this self-energy which produces the Xd/2 contribution. 

A more careful calculation for the regime d 3> A was performed by Matveev 
and Larkin , who considered the renormalization of A due to "pairing fluctuations" 
about the uncorrelated Fermi ground state \p/2)q. Summing up the leading logarithmic 
vertex corrections |Q [see inset of Fig. fH|(a)], they obtained a renormalized coupling 
A given, with logarithmic accuracy, by 



X + Xd 



u 
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1 - XdJ2 



1 
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2tt [iw - {sj - /i)] [-iw - (sj - /i)] 
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(45) 
(46) 
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This result evidently is valid only if d 3> uj^e^ 1 ^ ~ A/2 (which, incidentally, is 
another way of seeing that d ~ A defines the crossover between the fluctuation- 
dominated and bulk regimes). Matveev and Larkin concluded that 

(A^ L ) rcn ~ Ad/2 = d/ (2 log d/A) for d>A. (47) 

This logarithmic renormalization is beyond the reach of the g.c. variational BCS 
method, but was confirmed using exact methods |29| [42|, [35| (see Sec. |ll|). Its occur- 
rence, in a regime that in g.c. variational calculations appears to be "uncorrelated" , can 
be regarded as the "first sign of pairing correlations", in particular since, by Eq. (|46|), 
the interaction strength increases upon renormalization only if the interaction is at- 
tractive (A < would imply |A| < |A|). The pairing fluctuations responsible for this 
renormalization will be discussed in more detail in Sees, [ll] and |l2|. 

Unfortunately, A^ L is at present not measurable in ultrasmall gra ins, for the same 
experimental reasons as apply to £\/2 — £o (see footnote | in Sec. |9.l| ). 

9.4 Parity effect for pairbrcaking energies 

Braun and von Delft |3l], |1| Q discussed yet another parity effect, based on 

A P b = fi - n e , (48) 

the difference between the ■pair-breaking energies of an even and an odd grain [see 
Eq. (|42|)]. Fig. |l3](b) shows its behavior as function of d. In the bulk limit £l e ~ fl a ~ A 
and Ap b ~ 0. The most interesting feature of A p b is that it initially becomes negative 
as d/A increases; this occurs because in an odd grain pairing correlations are weaker 
and hence breaking a pair costs less energy than in an even grain. A p b becomes 
positive again for d/A > 0.5, since then pairing correlations are so weak that Ap b is 
governed by the kinetic energy cost of flipping a spin, which is 2c? for an odd grain but 
only d for an even grain. 

Ap b is directly measurable in RBT's grains, via the level-crossing fields /io,i = 
and = [Eq. ([4lj)]. The measured values fl e = 0.23 ± 0.01 meV and 

Q = 0.24 ± 0.01 meV cited in Sec. 3.1 give a positive value of Ap b = 0.1 meV, 



implying that the grain under study was too small to fall in the most interesting 
regime where A p b is negative. Braun and von Delft suggested that the latter should 
be observable in a somewhat larger grain with ft.1/2,3/2 < ho.i-, i-e. in Regime II of 
Fig. [□} (This suggestion assumes that despite the increased grain size, the com- 
plicating effect of orbital diamagnetism is still non-dominant in Regime II.) To look 
for negative Ap b experimentally would thus require good control of the ratio d/A, i.e. 
grain size. This might be achievable if recently-reported new fabrication methods, 
which allow systematic control of grain sizes by using colloidal chemistry techniques 
could be applied to Al grains. 
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Part II: Crossover from the bulk to the limit of a few electrons 

Part II of this review is devoted to the question: How do pairing correlations change 
when the size of a superconductor is decreased from the bulk to the limit of only a few 
electrons? In particular, we shall attempt to refine the answer given by Anderson [^jj , 
namely that superconductivity as we know it breaks down for d > A. 

First steps towards a more detailed answer were taken in the early 1970s by Stron- 
gin et al. |n| and by Muhlschlegel et al. (6^], who calculated the thermodynamic 
properties of ensembles small superconducting grains. Experimental realizations of 
such ensembles were, e.g., the granular films studied by Giaver and Zeller p0[ . 
The interest of theorists was rekindled in 1995 by RBT's success in probing individual 
superconducting grains. Apart from motivating the phenomenological theory of Braun 
et al. reviewed in part I, these experiments also inspired a substantial and still growing 
number of theoretical studies of how superconducting pairing correlations in 

such grains are affected by reducing the grains' size, or equivalently by increasing its 
mean level spacing d oc VoP 1 until it exceeds the bulk gap A. 

In the earliest of these, von Delft et al. studied the discrete BCS model of Sec. ^ 
within a parity-projected g.c. BCS approach [jl7| closely related to the variational BCS 
method of Sec. ^|. Their g.c. results suggested that pairing correlations, as measured 
by the pairing parameter or the condensation energy, vanish abruptly once d exceeds 
a critical level spacing dP^? that depends on the parity (p — or 1) of the number 

of electrons on the grain, being smaller for odd grains {dPjl? ~ 0.89A) than even 

grains (d^ cs ~ 3.6A). Parity effects were also found in a number of subsequent 
papers that used parity-projected g.c. methods to study the behavior of the BCS 
mean-field gap parameter A gc and related quantities as functions of level spacing |i~8|]- 
p8[, temperature and magnetic field. All these parity effects are consequences of the 



blocking effect (cf. Sec. 4.4): for odd grains, the unpaired electron somewhat disrupts 
the pairing correlations of the remaining paired ones, by reducing the phase space 
available for pair scattering. 

A series of more sophisticated canonical approaches [^9)-|37j (summarized in Sec. [Tl| ) 
confirmed the parity dependence of pairing correlations, but established that the abrupt 
vanishing of pairing correlations at dP^? is an artifact of g.c. treatments: pairing cor- 
relations do persist, in the form of so-called fluctuations, to arbitrarily large level spac- 
ings (24j, and the crossover between the bulk superconducting (SC) regime (c? <C A) 
and the fluctuation-dominated (FD) regime (d ^ A) is completely smooth pl-fSG]. 
Nevertheless, these two regimes are qualitatively very different |3l]]-|36|: the condensa- 
tion energy, e.g., is an extensive function of volume in the former and almost intensive 
in the latter, and pairing correlations are quite strongly localized around the Fermi 
energy ep, or more spread out in energy, respectively. Very recently, Di Lorenzo et 
al. |40) suggested that the remnant pairing correlations in the FD regime might be 
detectable via susceptibility measurements. 

Toward the end of 1998 and after the appearance of most of these works, R.W. 
Richardson pointed out (9^] to their various authors that the discrete BCS Hamiltonian 
on which they are based actually has an exact solution, discovered by him in 1963 Jl3"[ 
(and independently by Gaudin in 1968 |9q]). Richardson published his solution in the 
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context of nuclear physics in a series of papers between 1963 and 1977 Q]-[|5l| which 
seem to have completely escaped the attention of the condensed matter community. 
Very recently, the model was also shown to be integrable |)4| ^5). The revival of this 
remarkably simple exact solution after such a long and undeserved period of neglect is 
perhaps one of the most important consequences of RBT's experimental breakthrough: 
Richardson's solution allows the elucidation and illustration by exact means of many 
important conceptual ingredients of the standard BCS theory of superconductivity, 
such as the nature of pairing correlations, the importance of phase coherence, the 
validity of using a mean-field approximation and a grand-canonical formulation for 
bulk systems, and the limitations of the latter approaches for ultrasmall systems. 
Moreover, it allows the exact calculation of essentially all quantities of interest for 
ultrasmall grains. 

We shall therefore start part II by discussing the exact solution [Sec. [l(J. We then 
summarize the other canonical approaches somewhat more briefly than they perhaps 
would have deserved had an exact solution not existed, and compare their results to 
those of the exact solution [Sec. [fl]]. Next we analyze the qualitative differences be- 
tween the bulk and FD regimes [Sec. |l2), then discuss the case of randomly (as opposed 
to uniformly) spaced energy levels Sj [Sec. and finally discuss finite temperature 
parity effects [Sec. |l4). Throughout part II we set /i = 0, since canonical treatments 
make no reference to a chemical potential. 



10 Richardson's exact solution 



In this section we summarize some of the central results of Richardson's exact solution 
of the discrete BCS model. 



1 0.1 General eigenstates 



Consider N — 2n + b electrons, b of which are unpaired, as in Sec. 4.4. According to 
the general discussion there, the nontrivial aspect of solving the model is finding the 
eigenenergies £ n and corresponding eigenstates \^ n ) [Eq. (|l3|)1 of the pair Hamiltonian 
[Eq. (fu]), in which we set /i = below] 

u 

ij 

in the Hilbert space of all states containing exactly n pairs &t = ct + ct_ of electrons, 
where j runs over the set of all unblocked single-particle levels, U = I\B [I is the set 
of all interacting levels, B the set of all blocked levels]. In general, degenerate levels 
are allowed in J, but are to be distinguished by distinct j-labels, i.e. they have = Ej 
for i 7^ j . 

Richardson showed that the sought-after eigenstates (with normalization 
(\J r n |\[' n ) = 1) and eigenenergies have the general form 

n n U 

\^ n )=Mj{Bt\Q), £„ = 5>„ with flt = ^ _J_ . (50 ) 

v— 1 y—1 j 
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Here J\f is a normalization constant and the n parameters E v (y = 1, ...,n) are a 
solution of the set of n coupled algebraic equations 

1 U 1 " 2 

which are to be solved (numerically, see App. B2 of Ref. |16|j ) subject to the restrictions 
^ E v \i [i ^ v . Richardson originally derived this remarkably simple result by 
solving the Schrodinger equation for the wave-function . . . , j n ) of Eq. (|TT|). A 

simpler proof, also due to Richardson [^6|, may be found in Ref. |3(| and in App. B.l 
of Ref. its strategy is to verify that (Hu — £ n )\^n) — by simply commuting Hjj 
past the B\ operators in (|50|). 

Below we shall always assume the £j 's to be all distinct (the more general case that 
degeneracies are present is discussed by Gaudin |9q]). Then it can be shown explicitly 
[ p3[ that (i) the number of distinct solutions of Eq. (|5l|) is equal to the dimension of the 
n-pair Hilbert space defined on the set of unblocked levels U, namely ( N 7 ") 1 where Nu 
is the number of unblocked levels; and (ii) that the corresponding eigenstates ( |50| ) are 
mutually orthogonal to each other, thus forming an eigenbasis for this Hilbert space. 
This can easily be understood intuitively, since there exists a simple relation between 
the bare pair energies 2Ej and the solutions of Eqs. ([si]): as A is reduced to 0, it follows 
by inspection that each solution {E\, . . . , E n } reduces smoothly to a certain set of n 
bare pair energies, say {2Ej ± , . . . , 2£ Jji }; this particular solution may thus be labeled 
by the indices ji,...,j n , and the corresponding eigenstate ( |50|) written as \^ n ) = 
\jl,...j n ). By inspection, its A — > limit is the state \ji, . . . j n )u,o = Il"=i & ]j ): 
thus there is a one-to-one correspondence between the sets of all states {\ji, . . . ,j n )u} 
and {\jx, ■ ■ -jn)u,o}- -^ u ^ * ne hitter constitute a complete eigenbasis for the n-pair 
Hilbert space defined on the set of unblocked levels U, thus the former do too. 



10.2 Ground state 

For a given set of blocked levels -B, the lowest-lying of all states | ^ n , B) of the form (|l0|) , 
say \Sfr n ,B)G, is obtained by using that particular solution \jx,...j n ) I0r which the 
total "pair energy" £ n takes its lowest possible value. The lowest-lying of all eigenstates 
with n pairs, b blocked levels and total spin s — 6/2, say \n, s)q with energy £®(n), is 
that \ty n ,B)G for which the blocked levels in B all contain spin-up electrons and are 
all as close as possible to £f, the Fermi energy of the uncorrelated iV-electron Fermi 
sea \Fjsr). The E v for the ground state |n, s}q coincide at A = with the lowest n 
energies 2s j (j = 1, . . . , n), and smoothly evolve toward (initially) lower values when A 
is turned on, a fact that can be exploited during the numerical solution of Eq. ( f5l| ) . As 
A is increased further, some of the E^s become complex; however, they always occur 
in complex conjugate pairs, so that £ n remains real Q. For details, see Ref. Q and 
App. B.2 of Ref. []l6|| , where some algebraic transformations are introduced that render 
the equations less singular and hence simplify their numerical solution considerably. 
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10.3 General comments 

Since the exact solution provides us with wave functions, it is in principle straightfor- 
ward to calculate arbitrary correlation functions of the form (^jJbtfrt . . . 6^6 ./|v?n}) by 
simply commuting all b's to the right of all 6^'s. However, due to the hard-core boson 
commutation relations ([r]) of the 6's, the combinatorics is rather involved. Neverthe- 
less, Richardson succeeded to derive [EtJ explicit results for the normalization constant 
W of (|50| ) and the occupation probabilities u| and correlators Cjj of Eq. (|l^) (summa- 
rized in App. B.3 of Ref. [^6|). The exact result for the C^'s show that they are all 



positive, in agreement with the requirement (ii) formulated in Sec. 5.2. It is also natu- 
ral to ask whether in the bulk limit (d — > at fixed n d), the standard BCS results can 
be extracted from the exact solution. Indeed they can, as Richardson showed in |5l[] 
(following unpublished work by Gaudin |9j|), by interpreting the problem of solving 
the eigenvalue equations ( |5l| ) for the E u as a problem in two-dimensional electrostatics 
(see App. B.2 of Ref. |ll|). Exploiting this analogy, he showed that in the bulk limit, 
Eqs. (^lj) reduce to the well-known BCS equations determining the gap and chemical 
potential at T = [Eqs. ( |38| ) and (|39|)], and the ground state condensation energy 
£§ ond (n) [Eq. |[|)] to its BCS result, namely -A 2 /2d. 

Finally, let us mention that the Cambiaggio, Rivas and Saraceno have recently 
shown that the discrete BCS model is integrable and have constructed explicit ex- 
pressions for all its constants of the motion j94j. The latter's relation to Richardson's 
solution was clarified by Sierra j9j|, who has also explored possible connections be- 
tween the exact solution and conformal field theory. It would be an interesting chal- 
lenge for mathematical physicists to try to exploit this integrability to calculate finite- 
temperature properties exactly — although these can in principle be obtained from 
Richardson's solution by "simply" computing the partition function over all states, 
this is forbiddingly tedious in practice for large temperatures, since the eigenenergy of 
each state requires a separate (non-trivial) numerical calculation. 



11 Comparison of other canonical methods with the exact solution 

In this section we briefly mention the various canonical methods by which the discrete 
BCS model had been investigated prior to the revival of Richardson's exact solution 
in 1999. All of these studies used a half-filled band with fixed width 2cjd of uniformly- 
spaced levels [i.e. Ej = j d+ (1 — p)d/2, as in Eq. (0)], containing N = 2n + b electrons. 
Then the level spacing is d = 2u}j^/N and the bulk gap is A = £Jd/ sinh(l/A). Fol- 
lowing |3l[], we take A = 0.224 throughout this section. To judge the quality of the 
various approaches, we compare in Fig. |lj the results which they yield with those 
from Richardson's solution, for the even and odd (s = 0, 1/2) condensation energies 
£J™ nd _and the Matveev-Larkin parity parameter A^ L [cf. Sec. |9.3| . In the notation of 
Sec. 10.2| , these are given by 



ET nd (n) = £?(n)-(F N \fr\F N ), (52) 
A^(n) = £? /a (n)-[£ G (n)+£ G (n + l)]/2. (53) 

Following the initial g.c. studies [if, 111; @ 01 the discrete BCS model, the 
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first purely canonical study was that of Mastellone, Falci and Fazio (MFF) [£9|, who 
used Lanczos exact diagonalization. Despite being limited to n < 12, they managed 
to reach reasonably small ratios of d/A by using an ingenious scaling approach: for 
a given level spacing d, they increased the coupling constant A to about 0.5, thereby 
decreasing d/A = rf/oJD sinh(l/A) to values as small as 0.5. This allowed them to 
probe, coming from the few-electron side, a remarkably large part of the crossover to 
bulk limit. They found, La., that the condensation energies are negative for all d, 
showing that the system can always gain energy by allowing pairing correlations, even 
for arbitrarily large d. 

Berger and Halperin (BH) |50| showed that almost identical results can be achieved 
with less than 6 pairs, thus significantly reducing the calculational effort involved, by 
first performing a "poor man's scaling" renormalization: they reduce the bandwidth 
from o>d ~ nd to, say, wd ~ nd (with n < 6) and incorporate the effect of the removed 
levels by using a renormalized coupling constant, 



A = A 



^ 2e 

WD<|Ej|<"D 



(54) 



The reduced system is then diagonalized exactly. Note that the renormalization of 
Matveev and Larkin [Eq. (|46|)1 corresponds to taking cDd — d in Eq. (j54j), i.e. to 
integrating out the entire band. Also note that the renormalization prescription of 
( p4| ) has the property that it would leave the bulk gap invariant in the limit d/A — > 0, 
for which Eq. ( ]54| ) would imply uiue^ 1 ^ ~ ujue^ 1 ^ — A. 

To access larger values of n, Braun and von Delft fjlf used a fixed-n projected 
BCS approach (PBCS), in which BCS-like variational wavefunctions are projected to 
fixed particle number, as in Eq. (|30|). The projection integrals occurring in Eq. ( |30| ) 
were evaluated numerically for n < 600, using tricks developed in the nuclear physics 
literature by Bayman ^7|, Dietrich, Mang and Pradal |9£| and Ma and Rasmussen 
[ p9[ , and summarized in part in the book of Ring and Schuck |7(| . (A much simpler 
way of dealing with the projection, using recursion relations, was recently found by 
Dukelsky and Sierra j33|.) The PBCS method gives condensation energies that (i) 
are significantly lower than the grand-canonical ones [see Fig. thus the projection 
much improves the variational Ansatz, and that (ii) are negative for all d, confirming 
that the abrupt vanishing of the g.c. condensation energies is indeed an artifact of 
the g.c. treatment. The PBCS method is able to fully recover the bulk limit, but 
the crossover is not completely smooth and shows a remnant of the g.c. breakdown 
of pairing correlations: the d-dependence of the condensation energy (£;cond-jPBCS 
changes rather abruptly [kinks in the short-dashed lines in Fig. [l4|(a)1 from being 
extensive (~ to being practically intensive (almost d independent). 

It should be mentioned here that a generalization of the PBCS method to finite 
temperatures has been worked out by Essebag and Egido in the context of nuclear 
physics [jlOCfl . The PBCS method has recently also been applied to the attractive 
Hubbard model in one dimension by Tanaka and Marsiglio p^| , who found even-odd 
and super-even effects. The latter consist of differences between the number of pairs 
being equal to n = 2m or 2m + 1, and arise if boundary conditions are used that 
produce doubly-degenerate levels (eg = e_g) near the Fermi surface. 
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Fig. 14 (a) The even and odd (s = 0,1/2) condensation energies E c s ° nd of Eq. (||) [in 
units of A], calculated with BCS, PBCS and exact wave functions pa^ , as functions of 
d/A = 2sinh(l/A)/(2n + 2s), for A = 0.224. For comparison, the dotted line gives the 
"bulk" result Eo ulk = -A 2 /(2d). (b) Comparison || of the parity parameters Ap L @ of 
Eq. ( p3[ ) [in units of A] obtained by various authors: ML's analytical result (dotted lines) 
[A(l - d/2A) for d < A, and d/21og(ad/A) for d > A, with a = 1.35 adjusted to give 
asymptotic agreement with the exact result]; grand-canonical BCS approach (dash-dotted 
line) [the naive perturbative result ^Xd is continued to the origin]; PBCS approach (short- 
dashed line); Richardson's exact solution (solid line); exact diagonalization and scaling by 
MFF (open circles) and BH (long-dashed line). 



Dukelsky and Sierra (3^, |3^| used the density matrix renormalization group (DMRG) 
(with n < 400) to achieve significant improvements over the PBCS results for the dis- 
crete BCS model, in particular in the regime of the crossover, which they found to be 
completely smooth. In general, the DMRG approach is applicable to systems that can 
be divided into two pieces, called block and environment, which interact via a prefer- 
ably rather small number of states. One starts with a small block and environment, 
computes their combined density matrix, then enlarges both and recomputes the den- 
sity matrix, etc, until a large part of the system has been treated. Dukelsky and Sierra 
chose the block and environment to consist, respectively, of all particle or hole states 
relative to the Fermi sea (for a detailed description of the method, see |53|). Since the 
pairing correlations involving coherent superpositions of particle and hole states are 
peaked in a rather small regime of width A around the Fermi energy [compare Figs. |?j 
or n6[, the "interaction" between block and environment is "localized", so that the 
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DMRG can a priori be expected to work rather well for this problem. 

Finally, Dukelsky and Schuck |37j showed that a self-consistent RPA approach, 
which in principle can be extended to finite temperatures, describes the FD regime 
rather well (though not as well as the DMRG). 

To check the quality of the above methods, Braun ]42| |35[ | computed E^ ond (for 
s = 0, 1/2) and A^ L using Richardson's solution (Fig. |14[). The exact results 

(i) quantitatively agree, for d — > 0, with the leading — A 2 /2d behavior for E^ ond 
obtained in the g.c. BCS approach |l7|, [HI which in this sense is exact in 
the bulk limit, corrections being of order <?; 

(ii) confirm that the even ground state energy always lies below the odd one (this 
had independently been proven rigoroulsy by Tian and Tang J38|); 

(iii) confirm that a completely smooth []32| , |3^| crossover occurs around the scale 
d ~ A at which the g.c. BCS approach breaks down; 

(iv) show that the PBCS crossover |H] is qualitatively correct, but not quantitatively, 
being somewhat too abrupt; 

(v) are reproduced remarkably well by the approaches of MFF |2{| and BH (3(J ; 

(vi) are fully reproduced by the DMRG of (3^, with a relative error of < 10~ 4 for 
n < 400; our figures don't show DMRG curves, since they arc indistinghuishablc 



from the exact ones and are discussed in detail in 




The main conclusion we can draw from these comparisons is that the two ap- 
proaches based on renormalization group ideas work very well: the DMRG is essen- 
tially exact for this model, but the band-width rescaling method of BH also gives 
remarkably (though not quite as) good results with rather less effort. In contrast, 
the PBCS approach is rather unreliable in the crossover region. To study general- 
izations of the discrete BCS model, e.g. using state-dependent couplings of the form 
rfy~!- Ai,-&J&„-, the DMRG would thus be the method of choice. 

12 Qualitative differences between the bulk and the few-electron regimes 

Does the fact that the exact condensation energy E c s ond is always negative, even for 
arbitrarily large d/A, mean that the system stays "superconducting" even if it is ar- 
bitrarily small? The answer is certainly no, since in the fluctuation-dominated (FD) 
regime, the pairing correlations are qualitatively different than in the bulk, supercon- 
ducting regime. In this section we shall try to make this statement more precise by 
analyzing the qualitative differences between the two regimes, with regard to the A- 
and d-dependence of E^ ond , and the behavior of the occupation probabilities v?. 

Fig. |l5|(a) shows, on a log-log plot, the d-dependence of the even condensation en- 
ergy EQ° nd (d). Note that even on the log-log plot, the crossover of the exact i?g ond from 
the bulk to the FD regime is completely smooth. According to Sierra and Dukelsky 
HH , the exact result for EQ° nd (d) can be fitted very well to the form 

E c ond (d) = -A 2 /(2d) - Vo (ln2)Lu D X 2 + lQ (Ad/2uj B )log(2Lu B /d) , (55) 
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Fig. 15 Log-log plot |42j] of some of the curves of Fig. D_4|(a) for the even condensation energy 
Eo ond [in units of A], for A = 0.224; its asymptotic - A 2 / (2d) behavior for d/A -» is shown 
by the dotted line. 



where r]Q and 70 are constants of order unity [ p3[ . The first term is extensive (oc Vol) 
and dominates in the bulk limit; its standard heuristic interpretation [ |58| is that 
roughly A/d levels (those within A of ep) are strongly affected by pairing, with an 
average energy gain per level of —A/2. The second term, which is intensive and 
dominates in the FD limit, is equal (up to the numerical factor rj ) to the result 
from second-order perturbation theory Q, namely (Ad) 2 Yl7=i ^j- n+i (2 £ » — 2£j) _1 . 
This subleading term's d-independence (which was anticipated in |ElJ |78|) may be 
interpreted by arguing that in the FD regime, the number of levels that contribute 
significantly to i?Q° nd is no longer of order A/d: instead, fluctuations affect all n ~ 
2uiT>/d unblocked levels within cjd of ep (this is made more precise below), and each 
of these levels contributes an amount of order — (Xd) 2 /d (corresponding, in a way, to 
its selfenergy). Finally, the third term contains the small parameter A/cjd and thus 
represents a very small correction. 



The A- and volume-dependencies of EQ° nd in Eq. (55) strikingly illustrate the qual- 



itative differences between the bulk and FD regimes: in the bulk regime, dominated 
by the first term, i?Q 0nd is nonperturbative in A (since A ~ 2ujDer 1 / x ) and extensive, 
as expected for a strongly-correlated state; in constrast, in the FD regime, dominated 
by the second term, £'g ond is perturbative in A and practically intensive (up to the 
weak logd dependence of the third term). 

Perhaps the most vivid way of illustrating the qualitative difference between the 
bulk and FD regimes is to study properties of the ground state wavefunction. We shall 
consider here the correlators |3l| 

C](d) = {b]b j ){b j b]), (56) 

which measure the probability that a level can be "both occupied and empty", and 
vanish identically for states without pairing correlations. For the discrete BCS model 



Cj identically equals (bjbj) — {b\b-) 2 — v 2 — vj [by Eqs. (|16|) and ( |19| )] , which measures 



the fluctuations in the pair occupation number of level j, and it vanishes for any 
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-5 5 10 15 20 

(e.-8 F )/A 

Fig. 16 The occupation probabilities C 3 of Eq. @ for d/A = 0, 0.27, 1.09, 2.17 and 4.34 
[42]. In all three figures, the thick solid lines give the d — bulk BCS result, whereas circles 
and stars represent C^-values evaluated for discrete j's using the exact solution and PBCS 
method, respectively. For d = 0.27A, the PBCS and exact results are indistinguishable, and 
are shown in (a) as a single long-dashed line, which is also virtually identical to the bulk 
curve. For small d, pairing correlations are evidently localized within a few A of ep. With 
increasing d more and more weight is shifted away from £f into the tails; compared to the 
exact results, the PBCS method somewhat overemphasizes this derealization, which is one 
of the reasons why it produces a somewhat too abrupt crossover. 



blocked single-particle level. Note that also equals (bjbj) — {c*- + c- + ){^_c-_)\ 

this form, which was used in and corresponds to the diagonal terms under the sum 
in Eq. ( p2[) for A can , can be interpreted as the probability enhancement for finding a 
pair of electrons instead of two uncorrelated electrons in a single-particle level \j, ±). 

When evaluated using the grand-canonical BCS wavefunction, (C?)bcs is equal to 
UjVj = jA 2 /(e 2 + A 2 ) [thick solid lines in Fig. [l6], the same function as that plotted 
in Fig. 0] . The (C^bcs's thus have a characteristic peak of width oc A around Ep, 
implying that pairing correlations are "localized around £f in energy space" , which 
may be taken to be the defining property of "BCS-like correlations". Moreover, in 
the bulk regime d <C A, the (Cj )bcs are virtually identical to (C?j)exact [long-dashed 
line of Fig. |l^(a)], vividly illustrating why the grand- canonical BCS approximation is 
so successful: not performing a canonical projection hardly affects the parameters Uj 
and fry if d <C A, but tremendously simplifies their calculation. 

As one enters the FD regime d > A, the character of the correlator (Cj) exact 
changes [Fig. |l6|(b), circles]: weight is shifted into the tails far from ep a t the expense of 
the vicinity of the Fermi energy. Thus pairing correlations become delocalized in energy 
space (as also found in [ 
quite appropriate. In the extreme case d ^> A, the (Cj) eX act for all interacting levels 



29| 32l p3|), so that referring to them as mere "fluctuations" is 
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Fig. 17 Exact even and odd condensation energies E^° nd [in units of A] for equally spaced 
levels (dashed line), and the ensemble-average {Es° n ) for randomly-spaced levels (solid line) 
jj^ j. The height of the fluctuation bars gives the variances 5Eg° nd . 



are roughly equal. 

Richardson's solution can also be used to calculate, for a given set B of blocked lev- 
els, the rf-d epen dence of the canonical order parameter Af an (d) of Eq. (^). Schechter 
has found [ 1 1 1 that it can be fit to the form Af an (d) — A(l + 73d/ A), where 
is a positive numerical constant, and the linear term essentially reflects the factor of 
d in the definition of A^ n . The fact that A^ n is a strictly increasing function of d 
is in very striking contrast to the behavior of the grand-canonical pairing parameters 
A s (d) shown in Fig. |](a). 



13 Effect of level statistics 

Smith and Ambegaokar investigated the effect of level statistics on the crossover be- 
tween the bulk and FD regimes pc|| . In contrast to the uniform level spacing used 
in previous works, they employed a random spacing of levels, distributed according 
to the gaussian orthogonal ensemble. Using a g.c. mean-field BCS approach, they 
found, interestingly, that randomness enhances pairing correlations: compared to uni- 
form spacings (u.s.) jL7), it (i) on average lowers the condensation energy E^ ond to 
more negative values, (E^ ond ) < £^ ond (u.s.), but (ii) these still are parity dependent, 
(EQ° nd ) < (£^°2 d ). These results can readily be understood intuitively: pairing corre- 
lations become stronger the higher the density of levels around s-p , where pair- mixing 
costs the least energy. When determining the amount of pairing correlations for a 
set of randomly-spaced levels, fluctuations that increase the level density near e-p are 
thus weighted more than those that decrease it, so that randomness enhances pairing 
correlations. 

Although the g.c. mean-field treatment of Smith and Ambegaokar breaks down for 
mean level spacings much larger than A, just as was the case in Jl7 
conclusions (i) and (ii) are robust. Indeed, these were recently confirmed by Sierra 
et al. pfjfl, who used Richardson's exact solution to calculate E c ° nd for ensembles of 
random levels [Fig. [l7]]. Moreover, they found that the blocking effect responsible for 
(ii) manifests itself in the fluctuations too, which likewise are parity dependent: for 
example, Fig. [l?] shows that both the variances SE c s ond = [((E c s ond ) 2 } - (£™nd)2]i/2 



. 18, 19], their main 
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and the randomness-induced changes in condensation energies |(££ ond ) — -Eg° nd (u.c) | 
were larger for even than for odd grains. 



14 Finite temperature parity effects 

Although finite-temperature studies of the discrete BCS model are not of direct rele- 
vance for spectroscopic measurements of the BRT-typc (a finite T would simply smear 
out the discrete spectra, thereby blurring their most interesting features), they are 
important in their own right for extending our understanding of superconductivity in 
ultrasmall grains. Wc hence review several recent finite-T developments below. 

To begin, let us note that parity effects are of course not restricted to the T = 



limit discussed so far. To be observable p2[-|86 , they only require the temperature to 
be smaller than the free energy difference 8T ~ A — k&T\a{N c g(T)\ between an odd 
and even grain. Here N e g(T) is the effective number of states available for quasiparticle 



excitations at temperature T, and for d <C A is given by N c g(T) = y87rTA/ci 2 
[p2[. Below the corresponding crossover temperature where 6 T = 0, determined by 



k B T* r = Aln[7V cff (T c * )] and roughly equal to A/ In y8irA 2 /d 2 , the single unpaired 
electron begins to matter: it causes a crossover from e-periodicty to 2e-periodicity in 



the I-V characteristics of mesoscopic superconducing SET's |82| |83|, |84| |85|, |8J|, due 
to the ground state energy difference £i/ 2 — So — A. Since T c * becomes of order A in 
nanoscopic grains with d ~ A, parity effects should survive to temperatures as high 
as the (bulk) superconducting transition temperature T c itself. 

Regrettably, the canonical methods discussed in the preceding sections become im- 
practical at finite temperatures, since the number of states that need to be considered 
increases rapidly for T > d, A. On the other hand, g.c. finite-T methods, some of 
which we review below, are, in principle, inherently unreliable for d > A. This ap- 
plies in pa rticular to the simplest of these, parity-projected mean-field theory [^9[ jg7| 
(Sec. 14.1 ) and certain variational generalizations thereof 0, |22| (Sec. 14. 2j ): they 



yield the same sharp phase transition as function of temperature for finite systems 
as for bulk systems, whereas on general grounds no sharp transitions are possible in 
finite systems. The reason for this problem is that they neglect fluctuations in the 
order parameter, which become very important in the transition region. The sharp 
transition is smoothed out once fluctuations are included. A rather efficient way of 



doing this is the so-called static path approximation (Sec. 14.3). Its use is illustrated 



in Sec. 14.4 for a calculation of the spin susceptibility, which shows an interesting 



parity effect that should be measurable in ensembles of ultrasmall grains. 



14.1 Parity-projected mean-Bcld theory 

The simplest finite-T approach that is able to keep track of parity effects is parity- 
projected mean-field theory, first used in nuclear physics by Tanabe, Tanabe and 
Mang Jl02| , and, indepedently, introduced to the condensed-matter community by 
Janko, Smith and Ambegaokar [Q and Golubev and Zaikin J87j. One projects the 
g.c. partition function exactly onto a subspace of Fock space containing only even or 
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Fig. 18 d-and T-dependence of the pairing parameter A p / 2 (d,T), calculated using parity- 
projected mean-field theory |l7j]. Curve A gives the bulk gap A(0,T), with A(0,0) = A; 
curves B-E give A p / 2 (d, T)/A as a function of d/A and T/A for p = (B, D) and p = 1 
(C, E). The critical spacings do CS = 3. 56 A and d®f 2 s = 0.890A given here differ somewhat 
from those in Fig. ^(a) , because the present mean-field approach differs in minor details (via 
terms that vanish when d — > 0) from the variational approach of Sec. 3.2 



odd (p = 0, 1) numbers of particles, using the parity-projector P p : 

Zf = Tr KC P p e-^ 6 -^ , P ,i = |[1±(-1)*] . (57) 

One then makes the mean-field replacement bj — > {bj — (bj) p } + (bj) p , neglects terms 
quadratic in the fluctuations represented by { }, and diagonalizes H in terms of 
the Bogoljubov quasiparticle operators 7 J(J of Eq. (p8f). The self-consistency condi- 
tion Ap/2 = ^dJ2j(bj)pi evaluated in a parity-projected g.c. ensemble according to 
Eq. (p7p, leads to a gap equation of the standard form, 

J= d E ^(l-E^)' E j = y/( Ej -f l y + \A p/2 \^ (58) 

\e,\<ivo 3 V a / 

which is parity-dependent, via the occupation function f p j a — ('Jj a 'Yjcr)p f° r quasi- 
particles. Since their number parity is restricted to be p, f P j a differs from the usual 
Fermi function /? . The condition 2n+p = (N) p fixes the chemical potential \i to lie 
exactly half-way between the last filled and first empty levels if p = 0, and exactly on 
the singly-occupied level if p = 1, implying /i = in both cases [by Eq. (Q)]. 



von Delft et al. 17 1 applied this approach to the discrete BCS model with uniformly- 
spaced levels, and solved Eq. ( |58| ) for the parity parameter A p / 2 {d,T) as function of 
both level spacing and temperature. Fig. 18] summarizes their results. At zero tem- 
perature, A D / 2 ((i, 0) corresponds to the spin-dependent parity parameters A s=p j2 dis- 
cussed in Sec. |6.2| [cf. Fig. 0(a)], and drops to zero at a critical level spacing df™. The 
A p / 2 — > limit of Eq. ( |58| ) defines the parity-dependent "critical temperature" T c ^ p (d), 
which can be viewed as another measure of how rapidly pair-mixing correlations break 
down as function of level spacing (although ultrasmall grains of course cannot undergo 
a sharp thermodynamic phase transition, which can only occur if n — > oo). In both the 
even and odd cases, the behavior of T c p{d) shows direct traces of the parity projection: 
In the even case, T C) o(d) [Fig. [l^, curve D] is non-monotonic as function of increas- 
ing d, initially increasing slightly before dropping to zero very rapidly as d — > d^ cs . 
The intuitive reason for the initial increase is that the difference between the actual 
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and usual quasiparticle occupation functions is f P j a — f- a < for an even grain (be- 
coming significant when d ~ A), reflecting the fact that exciting quasiparticles two 
at a time is more difficult than one at a time. Therefore the quasiparticle-induced 
weakening of pairing correlations with increasing T will set in at slightly higher T if 
d~A. 

In the odd case, the critical level spacing d^£ s (T) [Fig. [l^, curve E] is non- 
monotonic as a function of increasing T, first increasing to a maximum before be- 
ginning to decrease toward d®^(T c ) — 0. The intuitive reason for this is that for 
< Aj/2 <C T,d, the odd j = function f p oa(T) becomes somewhat smaller than 
its T = value of 1/2, because with increasing T some of the probability for find- 
ing a quasiparticle in state j "leaks" from j = to higher states with j ^ 0, for 
which EJ 1 < Eq 1 in Eq. (pq). Thus, the blocking-of-pair-scattering effect of the odd 
quasiparticle becomes slightly less dramatic as T is increased, so that d^£ s increases 
slightly. 

It should be noted, however, that although the non-monotonicities of T C: o(d) and 
djS (T) are intuitively plausible within the g.c. framework in which they were derived, 
their physical significance is doubtful, since they fall in the regime where d/A s > 1 
and the g.c. approach is unreliable, due to its neglect of fluctuations. 



14.2 Variational extensions of BCS theory 

The above-mentioned results of von Delft et ol. Jl7| were reproduced and extended 
to finite magnetic fields by Balian, Flocard and Veneroni, using a more general g.c. 
variational BCS approach ED, |22j. It is designed to optimize the characteristic func- 
tion <p(£) = In TrP f) e~ /3 ^~' iiV ) A(g) , where P p is the parity projector of Eq. @ and 
j4(£) = exp(— J2 7 £-yQy), Q-y being observables of interest (e.g. the total spin) and 
^ 7 the associated sources (e.g. the magnetic field). This approach goes beyond the 
usual minimization of the free energy , since it optimizes not only thermodynamic 
quantities but also equilibrium correlation functions, which can be obtained by dif- 
ferentiating e^^' with respect to £ 7 . However, its g.c. version also suffers from the 
drawback of yielding abrupt, spurious phase transitions even though the systems are 
finite. Presumably this problem would be cured if an exact projection to fixed particle 
number were incorporated into this approach, but this is technically difficult and has 
not yet been worked out. 



14.3 Static path approximation 

For finite systems, in contrast to infinite ones, fluctuations of the order parameter 
about its mean-field value are very important in the critical regime, causing the phase 
transition to be smeared out; conversely, the spurious sharp transition found in the 
g.c. approaches above is a direct consequence of the neglect of such fluctuations. A 
rather successful way of including fluctuations is the so-called static path approximation 
(SPA), pioneered by Muhlschlegel, Scalapino and Denton |^2| and developed by various 
nuclear theorists [ 103 ]- [ 1 1 9 1 , while recently an exact parity projection has also been 
incorporated [Eq|-[E8|. A detailed and general discussion, including a complete list 
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Fig. 19 Temperature dependence of the pairing correlation energy A ca n of Eq. ( p9[ ) [in 
units of level spacing d], as calculated in [B6J using parity-projected mean-field BCS theory 
(dotted lines), SPA (dash-dotted), CSPA (dashed) and exact diagonalization (solid lines). 
A system of 10 equally-spaced, doubly-degenerate, levels was studied, containing 10 (left 
panel) or 11 (right panel) electrons. A can (T) is shown at four fixed values of d/A (thus this 
figure elaborates Fig. |li| ), namely 0.60, 0.91, 1.7, and 15, labeled by a,b,c and d, respectively 
(calculated using A = 0.55, 0.45, 0.35 and 0.2). CSPA data are shown only above the CPS A 
breakdown, which occurs at Tcspa < A/4 for the cases considered. The absence of dotted 
lines for the cases d (even) and c,d (odd) means that for these no nontrivial mean-field BCS 
solution exists. Of course, (A can ) cxact is nonzero nevertheless. The abrupt BCS transition 
is completely smeared out for the SPA, CSPA and exact results, for which the asymptotic 
decay at T > A can be shown H to be A can ~ |(A 3 d/T) 1/2 . 



of relevant references, was given very recently by Rossignoli, Canosa and Ring p6[ . 
We therefore confine ourselves below to stating the main strategies of the SPA and 
illustrating its capabilities by showing its results [Fig. for the quantity 

A c 2 an = (Ad) 2 ^[C y -(C y ) A=0 ]. (59) 

ij 

A can is reminiscent ^ of A can of Eq. ( [22] ) , and measures the increase in pairing correla- 
tion energy due to a nonzero coupling strength A. 

One starts by decoupling the quartic interaction H re & of Eq. (^) into a quadra- 
tic form using a Hubbard-Stratonovich transformation with a complex auxiliary field 
A(r) = A+(r) + iA~(r), with Matsubara-expansion A ± (r) = £ n A±e i2 ™ T /P in the 
interval r 6 (0, 0). The parity-projected partition function of Eq. ( |57| ) then has the 
following path integral representation (our notation is deliberately schematic; see p6[ ] 
for a precise version): 

Z p « fUdA+d&n 2 P [A], Z p [A]=Tr{ J P p re-^drMA(r)]| j (6Q) 



JThe definitions of Ac a n,aild A can differ by terms of order (rf/ uta ) 1 ' 2 ; for example, when evaluating 
both using |BCS) of Eq. jl7| ) and comparing to A gc of Eq. (|Lq), one finds (A can )BCS = A gc = 
(A can ) BC s +0[(d/u )D ) 1 / 2 ]. 
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j° 3 

The path integral can be treated at several levels of sophistication: 

(i) In the simplest, one uses a "fixed-phase saddle-point approximation" for the 
"static" n = modes and neglects all n ^ modes, i.e. one fixes the phase of Aq + 
IAq — |Ao|e 1( ^ by, say, setting cf) = 0, so that JdA^dAg is replaced by Jd|A |, 
and approximates this integral by its saddle-point value. The saddle-point condition 
for maximizing i? p [|Ao|] then yields the gap equation (^8|), thus this approach simply 



reproduces the parity-projected mean-field approach of Sec. 14.1, including its sharp 
phase transition (Fig. [ll| dotted lines). 

(ii) The next-best approximation is obtained if one writes J dA^dA^" = J Q 27r d0o 
J °° |Ao|d|Ao| and performs the phase integral fully. Remarkably, "liberating" the 
phase degree of freedom in this way already suffices to smooth out the phase transi- 
tion |117| , [26] , even if the / d|A | integral is again replaced by its saddle-point value, 
provided that the latter is found by now maximizing |Ao|2p[|Ao|] (i.e. including the 
factor |Ao| from the integration measure). This yields a modified gap equation with 
a nontrivial solution for arbitrarily large T, i.e. no abrupt transition. 

(iii) For finite systems, fluctuations about the saddle become large in critical re- 



gions. To obtain an improved descriptio n of t he la tter (Fig. 19, dash-dotted lines), the 



static path approximation (SPA) |62j, [ 103 1- [ 109 incorporates all static fluctuations 



exactly, via a (numerical) evaluation of the full integral J °° |Ao|d|Ao| d(f>o over all 
"static paths" . 

(iv) In the so-called correlated static path approximation (CSPA) (also called SPA- 
+RPA), small-amplitude quantum fluctuations around each static path are included 
too, by performing the remaining JdA^ Q integrals in the gaussian approximation 



[ 1 10 1- 116 , p5[|-[|28|. The CSPA yields qualitatively similar but quantitatively more 
reliable results (Fig. |l^, dashed lines) than the SPA, but breaks down below a tem- 
perature Xcspa, below which the fluctuations of the A n j modes become large at 
unstable values of |Ao|, causing the gaussian approximation to fail. 

(v) Finally, in the so-called canonical CSPA one projects the partition function 
not only to fixed number parity (as done throughout above) but also to fixed particle 
number, by performing a n int e grat i on over t he chemical potential (before performing 



any of the A^ integrals) [100, 107, 108, 109 1. However, this too is usually done only 



in the gaussian approximation (and would produce negligible corrections to the CPSA 
results for the quantities shown in Fig. |l9| ). 

Comparisons with exact diagonalization results fllgfl (Fig. [l9| solid lines) show that 
in its regime of validity (T > Tcspa), the CSPA produces results that are qualitatively 
completely similar and also quantitatively very close to the exact ones, whereas the 
quantitative agreement is significantly worse if only the SPA is used. Since the CSPA 
is conceptually simple, well-documented pt| and straightforward to implement, it 
seems to be the method of choice for not too low temperatures. A possible alternative 



is a quantum Monte Carlo evaluation of the path integral (60) [120, 121 1, but the 
numerics is much more demanding than for the CPSA, while the convergence at low 
T is in general rather poor, due to the familiar sign problem of Monte Carlo methods. 
The development of canonical finite-T methods that remain quantitatively reliable 
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for d > A and arbitrarily small T is one of the open challenges in this field. It 
would be very interesting if progress in this direction could be made by exploiting the 
integrability |94|, of the model, using Bethe Ansatz techniques. For the FD regime, 
another possibility would be to develop a finite- T generalization of the self-consistent 
RPA approach of Dukelsky and Schuck |57| . 



14.4 Re-entrant spin susceptibility 

For grains so small that (i > A, the spectroscopic transport measurements of BRT 
are not able, in principle, to reliably detect the effect of pairing correlations, since in 
this regime these cause only small changes to the eigenspectrum of a normal metallic 
grain, whose spectrum is, however, irregular to begin with. In contrast, thermody- 
namic quantities do have the potential to measurably reveal the existence of pairing 
correlations for d > A. Since very recently parity effects for the spin susceptibility 
have been observed experimentally for an ensemble of small, normal metallic grains 



122 , it is an interesting and experimentally relevant question to investigate how pair- 



ing correlations affect its behavior in superconducting grains. 



This question was worked out in detail by Di Lorenzo et at | 40 , whose results are 



summarized in Fig. |2Cj. The spin susceptibility for an isolated grain is defined as 
d 2 F p (T,H) 



X P (T) 



dH 2 



(62) 



H=0 



where T v = — k-QTXnZ^ 11 is the free energy of a grain with parity p and Zp an is the 
canonical partition function. 



In the bulk limit [Fig. 20(a)], it is well-known [123 1 that the spin susceptibility de- 
creases below its Pauli value XPauii = 2/zb 2 /d once T drops below the superconducting 
transition temperature T c , since the electrons tend to bind into Cooper pairs, which are 
spin singlets and do not contribute to the spin susceptibility. Interestingly, however, 
the spin susceptibility becomes parity-dependent as T is lowered below the crossover 



temperature T* r mentioned in the opening paragraphs of Sec. 14: In the even case, 
Xo exponentially drops to zero for sufficiently small temperatures, T <C max(A, d), 
for reasons that are intuitively obvious in the two limits A»d (all electrons bound 
into Cooper pairs) and A <C d (no Cooper pairs, but all levels doubly occupied). In 
contrast, in the odd case Xi shows a re-entrant behavior, in that it increases as ^b 2 /T 
for low temperatures, due to a Curie-like contribution from the unpaired odd electron. 
As a result, xi(T) nas a minimum somewhat below T* T , which can be viewed as a 
"smoking gun" for pairing correlations, since it is absent for odd normal grains. For 
these, Xi{T) [long-dashed A = curve in Fig. ^(b)] also has the Curie-like increase 
at very low T, but lacks the initial pairing-induced decrease as T is reduced below T c . 

Remarkably, Di Lorenzo et al. found that this reentrance of xi survives also for 
d > A [Fig. p0|(b)]: although pairing correlations survive here only as fluctuations, 
these are evidently sufficiently strong to still significantly reduce Xi{T) relative to 
XPauii [by several percent even for d/A ~ 50 (!)], before the Curie- like increase sets 
in at low T . Di Lorenzo et al. established this result by considering the limits T d 
and T <ti 2d analytically, using a static path approximation to capture the crossover 
numerically, and checking the results for T < d using Richardson's exact solution 
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Fig. 20 Spin susceptibility xo (xi) °f an even (odd) superconducting grain as function 
of T, plotted in units of its bulk high-T value XPauii = 2fiB 2 /d. (a) Bulk limit (d/A -C 
1): the even (solid) and odd (short-dashed) curves were calculated using parity-projected 
mean-field theory, the long-dashed curve using standard (unprojected) BCS theory, (b) 
Fluctuation-dominated regime (d/A > 1) for xi (main figure) and xo (inset): All curves 
were calculated using the static path approximation, except the short-dashed curves in the 
main figure, calculated using Richardson's exact solution, and the long-dashed curves for the 
non-interacting case (A = 0). 



(they considered all eigenstates with excitation energy up to a cutoff A ~ AOd, for 
grains with TV < 100 electrons). This check shows that the static path approximation 
somewhat underestimates the amount of pairing correlations (its minima for xi (T) are 
too shallow), but in general is in good qualitative agreement with the exact results, 
confirming that it is a useful and qualitatively reliable tool for describing the crossover 
regime. 



15 Summary and outlook 

The technique of single-electron-tunneling spectroscopy on ultrasmall metallic grains, 
applied to Al grains, has proved to be a very fruitful way of probing electron pairing 
correlations, and the way in which these are modified by level discreteness. It has, 
in particular, inspired theoretical attempts to quantitatively understand how pairing 
correlations change during the crossover from the bulk the limit of a few electrons. 
Let us briefly summarize the main conclusions reached in the preceding sections: 

Part I: For largish Al grains, the observation of a distinct spectral gap in even grains 
and its absence in odd grains is clear evidence for the presence of superconducting 
pairing correlations. These can be satisfactorily described using the simple discrete 
BCS model introduced in Sec. ||. The blocking of some levels by unpaired electrons 
leads to various measurable parity effects; among these, a pairbreaking-energy parity 
effect should be observable in experiments of the present kind, provided the grain size 
can be better controlled. In ultrasmall grains, the effect of a magnetic field on orbital 
motion is negligible. The dominant mechanism by which a magnetic field destroys 
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pairing correlations in ultrasmall grains is Pauli paramagnetism. Decreasing the grain 
size softens the first-order transition observed for thin films in a parallel field, by 
reducing the number of spins flipped from being macroscopically large for d <C A to 
being of order one for d ~ A. The grand-canonical variational BCS approach fails 
for d > A; nevertheless, it yields a useful framework for a qualitative analysis of the 
experiments, which had d < A. 

Part II: The crossover of the behavior of superconducting pairing correlations 
from the bulk limit (d <C A) to the fluctuation-dominated regime (d 3> A) is parity 
dependent and completely smooth. This remains true for systems with non-uniform 
rather than uniform level spacings. - Very remarkably, the discrete BCS model has 
an exact solution, due to Richardson, with which T = properties can be calculated 
rather easily. Finite-temperature properties for finite-sized systems can be calculated 
quite reliably with the correlated static-path approximation (provided T > Tcspa)- 
However, the development of canonical finite- T methods that remain quantitatively 
reliable for d > A is still an open problem. - The spin susceptibility x(T) of an odd 
grain shows an interesting reentrant behavior even for d 3> A, which might be a way 
to detect remnants of pairing correlations in the fluctuation-dominated regime. 

Finally, we would like to mention two further aspects of pairing correlations in 
ultrasmall metallic grains that have been discussed in the literature: the effect on the 
measured excitation spectrum of noncquilibrium excitations [j^, |^| , and of spin-orbit 
interactions |n[] . These topics have not been included here for lack of space, but have 
been reviewed in detail in sections 6.2 and 7.6 of Ref. Jlg| . 

Prospects for future work: 

Experiment: In the current generation of experiments, the grain's actual size and 
shape cannot be determined very accurately. It would be a great advance if fabrication 
techniques could be developed to the point that grains can be used which have been 
custom-made, by chemical techniques, to have well-defined sizes and shapes (e.g. spher- 
ical). This would significantly reduce the uncertainties which one presently encoun- 
ters when estimating characteristic parameters of the grain, such as the single-particle 
mean level spacing d or the Thouless energy. Moreover, it would allow systematic 
studies of the dependence of various quantities on grain size or mean level spacing [for 
example, it would be interesting to try to do this for the pairbreaking energies f2 e ,f2 
of Fig. |l^(b)]. Encouragingly, the feasibility of using chemically-prepared grains in 
SETs has already been demonstrated several times J89|, |9(], [n] , though the resulting 
devices have not yet been used for single-elcctron-tunncling spectroscopy. 

Theory: The behavior of superconducting pairing correlations in an individual 
ultrasmall grain can now be regarded as a subject that is well understood. It would 
be interesting to try to use the insights that have been gained for a single grain in 
order to now study systems of several coupled grains: what, for example, is the fate 
of the Josephson effect between two coupled grains as their sizes are reduced to the 
point that d ~ A? 
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